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Abstract 

We consider a class of nonlinear Klein-Gordon equations which are Hamiltonian and are 
perturbations of linear dispersive equations. The unperturbed dynamical system has a 
bound state, a spatially localized and time periodic solution. We show that, for generic 
nonlinear Hamiltonian perturbations, all small amplitude solutions decay to zero as time 
tends to infinity at an anomalously slow rate. In particular, spatially localized and time- 
periodic solutions of the linear problem are destroyed by generic nonlinear Hamiltonian 
perturbations via slow radiation of energy to infinity. These solutions can therefore be 
thought of as metastable states. The main mechanism is a nonlinear resonant interaction 
of bound states (eigenfunctions) and radiation (continuous spectral modes), leading to 
energy transfer from the discrete to continuum modes. This is in contrast to the KAM 
theory in which appropriate nonresonance conditions imply the persistence of invariant 
tori. A hypothesis ensuring that such a resonance takes place is a nonlinear analogue of 
the Fermi golden rule, arising in the theory of resonances in quantum mechanics. The 
techniques used involve: (i) a time-dependent method developed by the authors for the 
treatment of the quantum resonance problem and perturbations of embedded eigenvalues, 
(ii) a generalization of the Hamiltonian normal form appropriate for infinite dimensional 
dispersive systems and (iii) ideas from scattering theory. The arguments are quite general 
and we expect them to apply to a large class of systems which can be viewed as the 
interaction of finite dimensional and infinite dimensional dispersive dynamical systems, or 
as a system of particles coupled to a field. 
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1. Introduction and statement of main results - Smile of the Cheshire 

Cat* 

Structural stability and the persistence of coherent structures under perturbations of a dynami- 
cal system are fundamental issues in dynamical systems theory with implications in many fields 
of application. In the context of discrete, finite dimensional Hamiltonian systems, this issue 
is addressed by the celebrated KAM theorem |J, which guarantees the persistence of most in- 
variant tori of the unperturbed dynamics under small Hamiltonian perturbations. For infinite 
dimensional systems, defined by Hamiltonian partial differential equations (PDEs), KAM type 
methods have recently been used to obtain results on the persistence of periodic and quasiperi- 
odic solutions, in the case where solutions are defined on compact spatial domains with appro- 
priate boundary conditions 0] , [|HJ , . Variational methods have also been used to study this 
problem; see, for example, M. The compactness of the spatial domain ensures discreteness of 
the spectrum associated with the unperturbed dynamics. Therefore this situation is the gener- 
alization of the finite dimensional case to systems with an infinite number of discrete oscillators 
and frequencies. 

In this paper we consider these questions in the context of Hamiltonian systems for which the 
unperturbed dynamics has associated with it discrete and continuous spectrum. This situation 
arises in the study of Hamiltonian PDEs governing functions defined on unbounded spatial 
domains or, more generally, extended systems. The physical picture is that of a system which 
can be viewed as an interaction between one or more discrete oscillators and a field or continuous 
medium. In contrast to the KAM theory, where nonresonance implies persistence, we find 
here that resonant nonlinear interaction between discrete (bound state) modes and continuum 
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(dispersive radiation) modes leads to energy transfer horn the discrete to continuum modes. This 
mechanism is responsible for the eventual time-decay and nonpersistence of trapped states. The 
rate of time-decay, however, is very slow and hence such a trapped state can be thought of as a 
metastable state. 

The methods we develop are applicable to a large class of problems which can be viewed 
schematically in terms of "particle-field interactions". In the present work, we have not at- 
tempted to present general results under weak hypotheses but rather have endeavored to illus- 
trate, by way of example, this widely occurring phenomenon and clearly present the strategy 
for its analysis. A more general point of view will be taken up in future work. The approach 
we use was motivated and in part developed in the context of our study of a class of nonlinear 
Schrodinger equations with multiple nonlinear bound states and the quantum resonance prob- 
lem 1 60f1 , If)!]] , [62|, [ 63fl . See also [111,1121 and [fl6 |. Related problems are also considered in |32 



and 44 



We begin with a linear dispersive Hamiltonian PDE for a function u(x, t), x G !R n and t > 0. 
Suppose that this system has spatially localized and time-periodic solutions. Such solutions are 
often called bound states. A typical solution to such a linear system consists of (i) a non-decaying 
part, expressible as a linear combination of bound states, plus (ii) a part which decays to zero 
in suitable norms (dispersion). This paper is devoted to the study of the following questions: 

(1) Do small amplitude spatially localized and time-periodic solutions persist for typical non- 
linear and Hamiltonian perturbations? 

(2) What is the character of general small amplitude solutions to the perturbed dynamics? 

(3) How are the structures of the unperturbed dynamics manifested in the perturbed dynamics? 

Representative of the class of equations of interest is the nonlinear Klein-Gordon equation: 



d 2 t u = (A — V(x) — m 2 j u + Xf(u), AeR (1.1) 
with f(u) real- valued, smooth in a neighborhood of u = and having an expansion: 

f(u)=u 3 + 0(u 4 ). (1.2) 

Here, u : (x, t) i— > u(x,t) G 1R, x G 1R 3 , and t > 0. We shall restrict our large time asymptotic 
analysis to the case f(u) = u 3 . The more general case ( |1.2| ) can be treated by the technique of 
this paper by making suitable modification of W s,p norms used. 
We consider ( |1.1|) with Cauchy data: 

u(x, 0) = u (x), and d t u(x,0) = u\{x). (1.3) 
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Equation (|1 . 1| ) is a Hamiltonian system with energy: 

£[u,d t u] = \- J(d t u) 2 + |Vw| 2 + mV + V(x)u 2 dx + A J F(u) dx, (1.4) 

where F'{u) = f(u) and F(0) = 0. 

In the context of equations of type (1.1|), we have found the following answers to questions 
(1), (2) and (3). 

( Al) In a small open neighborhood of the origin, there are no periodic or quasiperiodic solutions; 
Corollary 1.1. 

(A2) All solutions in this neighborhood tend to zero (radiate) as t — > oo; Theorem 1.1. 

(A3) The time decay of solutions is anomalously slow*, i.e. a rate which is slower than the free 

dispersive rate; Theorem 1.1. 

Dynamical systems of the type we analyze appear in a number of physical settings. Consider 
a nonlinear medium in which waves can propagate. If the medium has local inhomogeneities, 
defects or impurities, these arise in the mathematical model as a spatially dependent coefficient 
in the equation (e.g. localized potential). Such perturbations of the original homogeneous 
(translation invariant) dynamics introduce new modes into the system (impurity modes) which 
can trap some of the energy and affect the time evolution of the system; see [|4|], |3"5| , |73|j . 



Let (K) = (1 + |fT| 2 ) 2 . For the nonlinear Klein-Gordon equation, (IT ), we prove the 
following result. 

Theorem 1.1. Let V(x) be real-valued and such that 
(VI) For 5 > 5 and \a\ < 2, \d°V(x)\ < C a (x)~ 5 . 

(V2) (-A + l)- 1 ((x • V) l V(x)) (-A + l)- 1 is bounded on L 2 for \l\ < N* with N* > 10. 
(V3) Zero is not a resonance of the operator — A + V; see J3l|/ , |72|/ . 
Assume the operator 

B 2 = -A + V(x) + m 2 (1.5) 

has continuous spectrum, a cont (H) = [m 2 , oo), and a unique strictly positive simple eigenvalue, 
Q 2 < m 2 with associated normalized eigenf unction, (p: 

B 2 ip = n 2 cp. (1.6) 

Correspondingly, the linear Klein-Gordon equation ( \1.1\ ), with A = 0, has a two-parameter 
family of spatially localized and time-periodic solutions of the form: 

u(x,t) = R cos(nt + 6) (p(x). (1.7) 
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Assume the resonance condition 



> 0. 



Here, P c denotes the projection onto the continuous spectral part of B and T c denotes the 
Fourier transform relative to the continuous spectral part of B. 

Assume that the initial data u , u\ are such that the norms \\u \\ W 2,2 nW 2,i and ||wi||wi,2 n wi,i 
are sufficiently small. Then, the solution of the initial value problem for (\1.1\), with A ^ and 



f(u) = u 3 decays as t — > ±oo. The solution u(x, t) has the following expansion as t — > ±oo: 

u(x,t) = R(t) cos (Qt + 9(t)) <p(x) + r)(x,t), (1.9) 

where 



3 , 



R(t) = 0(|t|-4) ,9(t) = <D(\t\*), and ||^(.,t)|| 8 = <D(\t\-*). (1.10) 
More precisely, 



R = R + 0(\R\ 2 ), ( \R\ small) 

t 

R(Q) = Rq, R 2 = \( V ,u )\ 2 + n- 2 \(<p, Ul )\ 



r = 2^r (l + ^rr 1 A 2 r \t\y*- (1 + o(\t\- 5 )) , 5>o 



Corollary 1.1. Under the hypotheses of Theorem 1.1, there are no periodic or quasiperiodic 
orbits of the How t > (u(t), d t u(t)) defined by ( \1.1\ ) in a sufficiently small neighborhood of the 
origin in the space (W 2 ' 1 n W 2 ' 2 ) x (W 1 ' 2 n W 1 ' 1 ). 

It is interesting to contrast our results with those known for Hamiltonian partial differential 
equations for a function u(x,t), where x varies over a compact spatial domain, e.g. periodic or 
Dirichlet boundary conditions [||], |T7|], | 57[| . For nonlinear wave equations of the form, ( |1 . 1| ) , 



with periodic boundary conditions in x, KAM type results have been proved; invariant tori, 
associated with a nonresonance condition persist under small perturbations. The nonresonance 
hypotheses of such results fail in the current context , a consequence of the continuous spectrum 
associated with unbounded spatial domains. In our situation, non-vanishing resonant coupling 
(condition ( |1.8| )) provides the mechanism for the radiative decay and therefore nonpersistence 
of localized periodic solutions. 
Remarks: 
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1) The condition (|1.8|) is a nonlinear variant of the Fermi golden rule arising in quantum 



mechanics; see, for example, [18|, [62]. This condition holds generically in the space of potentials 



satisfying the hypothesis of the theorem. Note that the condition (|1.8| ) implies that 

3ft G a cont {B) (1.11) 

i.e. the frequency, 3ft, generated by the cubic nonlinearity lies in the continuous spectrum of 
B. The hypothesis Q1.8Q ensures a nonvanishing coupling to the continuous spectrum. 

(2) The regularity and decay hypotheses on V(x) are related to the techniques we use to obtain 
suitable decay estimates on the linear evolution in L 2 ((x)^ cr dx) and LP] see section 2. 

(3) Persistence of dynamically stable bound states for special perturbations: There are nonlinear 
perturbations of the linear Klein-Gordon equation, with A = 0, for which there is a persis- 
tence of time-periodic and spatially localized solutions. Suppose we extend our considerations 
to the class of complex- valued solutions: u : IR™ x 1R — > C, and study the perturbed equation: 

(d? - A + V{x) + AImI^ 1 )^ = 0, (1.12) 



where 1 < p < oo for n — 1,2 and 1 < p < for n > 3. Unlike ( |1.1| ), equation ( |1.12j ) 



has the symmetry u \— > ue n and it has been shown []5T| that ( |1.12| ) has time periodic and 



spatially localized solutions of the form e lujt §(x; u) with $ G H x , which bifurcate from the 
zero solution at the point eigenvalue of —A + V(x) — uj 2 , by global variational [|39|, |65| and 



local bifurcation methods |45|]. There are numerous other examples of equations for which the 
persistence of coherent structures under perturbations is linked to the perturbation respecting 
a certain symmetry of the unperturbed problem. The stability of small amplitude bifurcating 
states of (1L12D can be proved using the methods of |24|], [fT0|l . An asymptotic stability and 



scattering theory has been developed for nonlinear Schrodinger dynamics (NLS) in [pQ] . If the 
potential V(x) supports more than one bound state, the above methods can be used to show 



the persistence of nonlinear excited states. However, it is shown in |63| that the NLS excited 
states are unstable, due to a resonant mechanism of the type studied here for (|1 . 1| ) . 
(4) In contrast with the time-decay rates associated with linear dispersive equations, the time- 
decay rates of typical solutions described by Theorem 1.1 is anomalously slow. See section 8 for 
a discussion of this point. 

We now present an outline of the ideas behind our analysis. For A = 0, solutions of equation 
( |1 . 1| ) are naturally decomposed into their discrete and continuous spectral components: 

u(x,t) = a(t)tp(x) + r}(x,t), ((p,r)(-,t)) = 0, (1.13) 
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where (/, g) denotes the usual complex inner product on L 2 . The functions a(t) and r](x,t) 
satisfy system of decoupled equations: 



with initial data 



a" + Q 2 a = 0, (1.14) 
d 2 r] + B 2 r] = 0, (1.15) 



a(0) = {<p,Uo), a'(0) = (y?,ui) 
r/(x, 0) = P c m , d t r)(x,0) = P c mi (1.16) 



For A 0, and for small amplitude initial conditions, we use the same decomposition, ( 1.13Q 



Now the discrete and the continuum modes are coupled and the dynamics are qualitatively 
described by the following model system: 

a" + Q 2 a = 3Xa 2 { X ,v{-,t)) (1-17) 

d 2 r]- Ar] + m 2 r] = X a 3 x- (1-18) 



Here, x( x ) is a localized function of x\ in particular, \ = (p 3 , and we assume (see ( |1 . 8| ) ) 3Q > m. 
In selecting the model problem ( |1.17| - |lTT8l) , we have replaced B 2 , restricted to its continuous 
spectral part, by B 2 = —A + m 2 , which intuitively should lead to the same qualtitative result. 

The system ( |1.17| - |iTT8| ) can be interpreted as a system governing the dynamics of discrete 
oscillator, with amplitude a(t) and natural frequency Q, coupled to a continuous medium in 
which waves, of amplitude r)(x,t), propagate, or as an oscillating particle coupled to a field. [] 

The system ( |1 . 1 7| - |iTT8"| ) is a Hamiltonian system with conserved total energy functional: 

1 f 1 

S[7],d t 7],a,a' } = - (d t r]) 2 + \Vr]\ 2 + m 2 r] 2 dx + - (a' 2 + tt 2 a 2 ) 

- Xa 3 / x(x)r)(x,t) dx (1-19) 



We now solve ( |1.18| ) and substitute the result into (|1.17j) . Note that, to leading order 



solutions of ( [1 . 1 7j ) oscillate with frequency Q. Therefore the a 3 term in ( 1.18Q acts as an external 

driving force with a 3Q frequency component. Since 9fl 2 is larger than m 2 , a nonlinear resonance 

of the oscillator with the continuum takes place. To calculate the effect of this resonance requires 

a careful analysis involving (i) a study of singular limits of resolvents as an eigenvalue parameter 

1 A related example is a model introduced in 1900 by H. Lamb |3S|| , governing the oscillations of mass-spring- 
string system. See also §, ||, ||, |o). 
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approaches the continuous spectrum (see section 4) and (ii) and the derivation of a normal form 
which is natural for an infinite dimensional conservative system with dispersion (see section 5). 
This leads to an equation of the following type for a(t) (or rather some near-identity transform 
of it): 

a" + (tt 2 + 0{\a\ 2 ))a = -V r 4 a', t>0 (1.20) 

where r = 0(|a|), and T > is the positive number given in (|1.8|) . For the model system (|1.18| ), 
T c in ( |1.8|) is replaced by the usual Fourier transform. This is the equation of a nonlinearly 
damped oscillator: 

j ((a') 2 + [n 2 + 0(\a\ 2 )}a 2 ) = -2Tr\a!) 2 < 0, for t > (1.21) 

Therefore, nonlinear resonance is responsible for internal damping in the system; energy is lost 
or rather transferred from the discrete oscillator into the field or continuous medium, where it 
is propagated to infinity as dispersive waves. Solutions of ( |1.20[ ) decay with a rate C(t _1//4 ), 
and it follows that \\r](-, £)||oo = 0(t~ 3 ^). Note, however that the total energy of the system, an 
H l type quantity, is conserved. In ( |1.20|) , we have neglected higher order effects coming from 
the continuous spectral part of H. These, it turns out, has a small effect and can be treated 
perturbatively. 

Aspects of the analysis are related to our recent treatment of the quantum resonance problem 
PHI , |62|] . The situation with quantum resonances can be summarized briefly as follows. Let Hq 
be a self-adjoint operator having an eigenvalue, A , embedded in its continuous spectrum, with 
corresponding normalized eigenfunction ip . Let W be a small localized symmetric perturbation, 
satisfying the (generically valid) Fermi golden rule resonance condition: 

r = vr (WiPo, 6(H - A )P«Wo) ± 0. (1.22) 

Then in a neighborhood of Ao, we show that the spectrum of H = H + W is absolutely 
continuous by proving that all solutions of the Schrodinger equation: 

id t i> = Htp = (H + W)ij, (1.23) 

with initial data which is spectrally localized (with respect to H) in a neighborhood of Ao, decay 
to zero in a local energy norm as t — > ±oo. 

Such solutions are characterized by exponential time-decay for an initial transient period, and 
then algebraic (dispersive decay) thereafter. During this transient period, A(t), the projection 
onto the mode ipo, is governed by the equation: 

A\t) = (zA* - r )A(t), t > (1.24) 
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where A* ~ A + (^o, Wip Q ), and by flOg) , T > 0. 

In the class of nonlinear problems under consideration, if we express the amplitude, a(t), as 



a(t) = A(t) e lUt + A(t) e~ im (1.25) 
then we find after a near-identity transformation an equation of the form: 

A' = ic 2 i\A\ 2 A + (zc 32 -^r)|A| 4 A t>0. (1.26) 



From this, the 0(t~ 1/4 ) behavior is evident. Equations ( |Q5| - |lT2^ ) lead to ([L20| ). As in (|T20|) , 



we have in ( |1.24|) and ( |1.26|) neglected the higher order coupling to the continuous spectral 



(radiation) component of the solution. These contributions are treated perturbatively. 
Remarks: 

(1) Lamb shift: Note that in the nonlinear problem, asymptotically there is no Lamb shift type 
correction to the frequency; the frequency shift is (9(|t4| 2 ) = 0(t~ 1 / 2 ) as t — > ±oo. 

(2) Emergence of irreversible behavior from reversible dynamics; dissipation through dispersion: 
Being Hamiltonian, the underlying equation of motion, ( [Lip , is time reversible. In particular, 
the equation has the invariance: u(x,t) i— > u(x, —t). Yet, the equation in ( |1.26| ) is clearly not 



time-reversible. This apparent paradox is related to the " ' e- prescription" discussed in section 4 
and Proposition 2.1; the singular limit 



lim exp(zV-A + m 2 t)(-A + m 2 - E ± ie)' 1 (1.27) 

satisfies a local decay estimate as t — »■ ±oo. For t < 0, the corresponding equation for A(t) 
would have -r replaced by +r ; cf. @, 0, 0, @. 



The nonpersistence of small amplitude spatially localized and time-periodic or quasiperiodic 
solutions to ( |1.1| ) has been studied in . p71 . In this work, the phenomenon of nonpersistence 
is formulated as a question concerning the instability of embedded eigenvalues of a suitable 
linear self-adjoint operator. A result on structural instability is proved; under the hypothesis 
( |1.8|) , time-periodic or quasiperiodic solutions of the linear problem equation (A = 0) do not 
continue to nearby solutions of the nonlinear problem (A ^ 0). 

The question of nonpersistence of periodic solutions has also been considered extensively in 
the context of the sine-Gordon equation 

d 2 u = dlu-sinu. (1.28) 
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Spatially localized and time-periodic solutions of the sine-Gordon equation are called breathers 
and the question of their persistence under small Hamiltonian perturbations in the dynamics 
has been the subject of extensive investigations. See, for example, f54}] , fl6|l , [ p~3f| |69| , ||, [f|, 

H, gg, h, gg. 



Analytical, formal asymptotic and numerical studies strongly suggest that for typical Hamil- 
tonian perturbations of the sine-Gordon equation, for example, the 4 model: 



dlu -u + u 3 



;i.29) 



no small amplitude breathers exist and that solutions obtained from spatially localized initial 
data radiate to zero very slowly as t tends to infinity. We believe this is related to the 
mechanism for slow radiative decay, as explained by Theorem 1.1 in the context of (|1 . 1|) . 

Finally, we wish to comment on the connection between our work and the approach taken in 
||56|| , ||, fL9| , and |33 |. As in the continuation theory of periodic solutions of ordinary differential 
equations Jl5|, it is natural to seek a periodic solution of (|1 . 1|) for A ^ which behaves, as the 
amplitude a tends to zero, like a solution ( |1 . 7| ) of the linear limit problem. The equation is 
autonomous with respect to time, so we seek a 27rf2~ 1 -periodic solution for A ^ with an 
amplitude dependent period. Since we do not know the period a priori, it is convenient to 
define 

u(x, t) = U a (x,s), s = Q a t (1.30) 
and require that U a (x, s) be 2n periodic in s. Thus (|1.1| ) becomes 



(Ql d 2 + B 2 )u a = A Ul 



[121) 



We formally expand the solution and frequency: 



JJ a = allx + a 3 U 3 + 
h 



Q. a = + a 2 fi 2 + 



;i.32) 



Substitution of (|1.32j ) into ( |1.31| ) and assembling terms according to their order in a, one gets a 
hierarchy of equations beginning with 



0{a l ) 
0(a 3 ) 



Q 2 d 2 s + B 2 



n 2 d 2 



U x = 

u 3 = xuf - 2fifi 2 a 2 t/i 



;i.33) 
;i.34) 



2 In breather type solutions have been constructed for the discrete sine-Gordon equation, where is 
replace by its discretization on a sufficiently coarse lattice. Also, a generalization of the notion of breather has 
been considered in the geometric context of wave maps pa). 



10 



Equation (|1.33|) has a solution 

Ui(x,s) = cos s (p(x). (1.35) 



Substitution into ( 1.34 ) gives the following explicit equation for U 3 : 



(tt 2 d 2 s + B 2 ) U 3 = (—ip 3 + 2tttt 2 (p J coss +-<^ 3 cos3s (1.36) 



We now express U$ in the form £7 3 = 




(fl 2 d 2 + B 2 ) UP = + 2Qn 2 (p coss (1.37) 



(tfd 2 + B 2 ) E/f = -p 3 cos 3s. (1.38) 

Since the inhomogeneous term in ( |1.37| ) is nonresonant with the continuous spectrum of B 2 , one 
can find a 27r-periodic solution provided the right hand side is L 2 (S\^ x 1R 3 ) - orthogonal to 
the adjoint zero mode: coss tp(x). This latter condition uniquely determines Q2- Note however 
that the right hand side of ( |1.38|) is resonant with the continuous spectrum of B 2 if 3fl > m. To 
compute the obstruction to solvability, seek a solution of (|1.38| ) of the form = cos 3s F. 
Then, 

( B 2 - 9n 2 ) F = ^ip 3 (1.39) 

and thus we expect to find F G L 2 if and only if the component of ip 3 in the direction of the 
generalized eigenfunction at frequency 3f2 of the operator B vanishes. Therefore, if the nonlinear 
Fermi golden rule ( |1.8[ ) holds, our formal expansion in L 2 (S\^ x 1R 3 ) breaks down. The relation 
with our work is that we prove that this obstruction to solvability, in fact, implies the radiative 
behavior of solutions described in our main theorem. 
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2. Linear Analysis 



In this section we summarize the tools of linear analysis employed in this paper. 
2.1. Estimates for the (free) Klein Gordon equation 

We begin by considering the Cauchy problem for the linear Klein-Gordon equation for a function 
u(x,t), x e JR n , t^O. 

d 2 t u = (A - m 2 )u = -Blu (2.1) 

u(x, 0) = u (x), d t u{x, 0) = ui(x) . (2.2) 

There exist operators EQ°\t) and E[°\t) such that 

u(x,t) = EP(t) u + E[°\t) u t . (2.3) 

We write 

E$\t)f = cos5ot/,and 

E[°\t)g = ^g, (2.4) 



where (w(k) = \/m 2 + k 2 ): 

E^\t)f = J cos uj{k)t f(k) e ik ' x dk 
<0) _ f sm u(k)t ^ Jk . x 



E[ U) (t)g = I " ' g(k)e lk - x dk. (2.5) 



The first result we cite is proved using stationary phase methods; see ||. See |36[] for another 
approach to decay estimates for (|2.1| ). Let W s,p (JR n ) denote the Sobolev space of functions with 
derivatives of order < s in LP. 



Theorem 2.1. Let 1< p < 2 1 + A- = I, 8 = ± - ± and < 9 < 1. Then, 

* — ' p p 2 %>' — — ' 

(a) if 5(n + 1 + 8) < v + s, we have for u — 0,1: 

II E^g\\ pl < K(t) \\g\\ S;P ,t>0, where (2.6) 

K(t) = cr^ n - x - e)5 , o < t < i, 

= cr (n - 1+e)s , t > 1. (2.7) 

(b) If s > (1/2 — l/p')(n + 2) — 1, then fort > 1 the (Schrddinger-like) L p ' decay rate, t _n ^~p 7 ^ 
holds in @). 
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In subsequent sections, we shall use some specific consequences of this result. 
Corollary 2.1. Consider the linear Klein Gordon equation (271) in dimension n = 3. Then, 



\\E?\t)g\\ 8 < Ct-l |M| lif t> 1 (2.8) 
\\E{°\t)g\\ 8 < Crt \\ g \\ h s 0<t< 1 (2.9) 
||££ 0) (*)<7||4 < IMIi,§, *>0, v = 0,1. (2.10) 

proof: Estimates ( |2.8[ ) and (|2.9| ) follow from the theorem with the choice of parameters: p' = 8, 
s = 1, n — 3, and = 1. Estimate (|2.10|) follows from the theorem with the choice of parameters: 



p' = 4, s = 1, n = 3 and = 0. 

2.2. Estimates for the Klein Gordon equation with a potential 

We now consider the Cauchy problem for the linear Klein-Gordon equation with a potential 

d 2 t u = (A - m 2 - V(x))u = -B 2 u (2.11) 

where B 2 is positive and self-adjoint. 

We write the solution to the initial value problem for Q2.ll ) with initial data ( |2.2| ) as: 



u(x,t) = E Q (t)u + Ei(t)ui , where (2.12) 

E (t)f = cos(Bt)f , and (2.13) 

^ / x sin(Bt) 
Ei(t)g = ~^9- 

One expects that estimates of the type appearing in Theorem 2.1 and Corollary 2.1 will hold 
as well for E (t) and Ei(t) restricted to P C L 2 , the continuous spectral part of H = B 2 . One can 
relate functions of the operator H, on its continuous spectral part, to functions of H = B 2 
using wave operators. Let 

W+ = strong - lim e itH e~ itHo . (2.14) 

t — >oo 

Wave operators relate Hq to H on P C L 2 via the intertwining property: 

H =W + H Wl on P C L 2 . (2.15) 
K. Yajima [72| has proved the W ,p (TR, n ) boundedness of wave operators for the Schrodinger 



equation. A consequence of this work is the following result for spatial dimension n = 3: 
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Theorem 2.2. Let H = —A + V, where V(x) be real-valued function on IR 3 satisfying the 
following hypotheses, which are satisfied by smooth and sufficiently rapidly decaying potentials: 
For any \a\ < £ there is a constant C a such that 



dx a 



(x) 



< C a {x) 



-6 



5 > 5 . 



Additionally, assume that is neither an eigenvalue nor a resonance of H; see jT^j, pi /. If 
< k, k! < I and 1 < p,p' < oo, then there exists a constant C > such that for any Borel 
function f on R 1 we have 

C 1 \\f(Ho)\\B(W k 'P,W k ''P') — \\f(.H)^ > c{H)\\B(W k <P,W k '' p ') — C \\f(Ho)\\ B ( W k,p tW k',piy 

Here, P C (H) denotes the projection onto the continuous spectral part of the operator H . 

Remark: In our applications, we shall use Theorem 2.2 with \l\ < 2. 

Theorem 2.2 implies that the dispersive estimates for the free Klein Gordon group carry over 
to the operators -E'o(^) an d Eiif) restricted to the continuous spectral part of B. 

Theorem 2.3. If g G Range P C (H), then each of the estimates of Theorem 2.1 and Corollary 
2. 1 hold with Ef^ replaced by E~ . 



2.3. Singular resolvents and time decay 

As discussed in the introduction, the damping term in the effective nonlinear oscillator (|1.26|) , 
is related to the evaluation of a singular limit of the resolvent as the eigenvalue parameter 
approaches the continuous spectrum. To ensure that the correction terms to the nonlinear 
oscillator (|1.26| ) can be treated perturbatively, we require local decay estimates for the operator 
e lBt (B — A ± z0) _1 , where A is a point in the interior of the continuous spectrum of B (A > m). 
Such estimates are analogous to those used by the authors in recent work on a time dependent 
approach to the quantum resonance problem [|6^j . 

We begin with a proposition, which is essentially a restatement of Theorem 2.3 for Ej(t), 
and is proved the same way. 

Proposition 2.1. (W s ' p estimates for e lBt ) Assume that V(x) satisfies the hypotheses of The- 
orem 2.2 (n = 3 ); see 172]] for hypotheses for general space dimension, n > 3. Let p and p' be as 
in Theorem 2.1 and let I > s > (| — -y)(n + 2), where I is as in Theorem 2.2. Then, 

II e lBt P c < C\t\- n{ ^ ] HVIU t ^ 0, (2.16) 
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The next proposition is the "singular" resolvent estimate we shall require in section 7. Let 
a*(n)=max{f,2 + ^}. 

Proposition 2.2. Assume the hypotheses of Proposition 2.1. Additionally, assume that V(x) 
satisfies hypothesis (V2) of Theorem 1.1. Also, let a > cr*(n). Then, for any point A (A > m) 
in the continuous spectrum of B we have for I = 0, 1, 2: 

\\{x)-° e iBt (B-A + zOy 1 P c (x)-° Vila < C (t)-& ||V||i,2, t > 0, 

|| e iBt (B - A — i0)' 1 P c (x)">|| 2 < C (t)-&* HVIIi.2, * < 0, (2.17) 

proof. We prove the estimate Q2.17| ) for the case of the £ > 0. For the case, t < 0, the same 
argument with simple modifications applies. 

Let g A = 9a{B) denote a smooth characteristic functions of an open interval, A which 
contains A and is contained in the continuous spectrum of B. Also, let ~g A = 1 — g A . We then 
write, for I — 1, 2: 

e *B* (S-A + ie) _I P c = e iBt g A {B) {B-A + ie)~ l P c + e iBt g A (B) {B-A + ie)~ l P c 

= sm + sm- (2.i8) 

We now estimate the operators Sj = lim^o Sp J = 1,2 individually. 

Estimation of Si (t) : 
Consider the case I = 2; the case I — 1 is simpler. First, we note that: 

(*>-* sm ( X )-° 

poo poo 

= -e it( - A ~ i£) / ds dr (x)- a e iT{B - A+l£) g A (B)P c (x)~ a (2.19) 

Jt Js 

Now we claim that the operator in the integrand satisfies the estimate: 



\\(x)- a e^ B - A+ ^ g A (B)P c (x)-1 B(L2) < C |— , (2.20) 

for r < a. The desired estimate on S[ then follows by use of ( |2.20| ) in ( |2.19|) and that a > o-*(n). 

It is simple to see that (|2.20| ) is expected. For, suppose that instead of BP C we had B . Let 
u(k) = \Jm 2 + k 2 , the dispersion relation of B . Then, setting L = (idk^u^j d^, and using 
that |VfcO;| 7^ on the support of g A , we get 

e iBot g A {B )f = J e iuj{k)t e lk ' x g A {k) j \k) dk 

= r r J U (V w{fc) ') e ik ' x g A {k)f{k) dk 

= t~ r f e luj{k)t (L^ r (e ikx g A {k)f{k)) dk. (2.21) 
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Estimation of (x) u e lBot g/\(Bo)f in L 2 by Fourier transform methods, using that | Vu>\ is bounded 
away from zero on the support of g&, we have that if a > max{|, 2} and r < a, then 



l<*> 



-cr iBot 



9A(B )f\\ 2 < C{t)- r \\{xYf\\ 



2- 



(2.22) 



Use of this estimate, in the above expression for S{ (with B replaced by Bq) gives, for I = 0, 1, 2: 

|| (x)-°e iBot (B - A + i0)~ l (x)- CT || B(L2) < C7<t>- r+I . (2.23) 



Two approaches to proving an estimate of the type flTgp and ( ^23|) for e iBt g A (B)P c are as 
follows: 

(a) One can "map" the estimate for S to that for BP C using the wave operator W + . In this 
approach we would need to derive estimates for W + g&(B ) in weighted L 2 spaces, or alternatively 

(b) One can use the approach based on the "Mourre estimate", developed in quantum scattering 
theory. This approach is more in the spirit of energy estimates for partial differential equations 
and does not require the use of wave operators. 

Here, we shall follow the latter approach. Time decay estimates like (|2.22| ) are a consequence 



of an approach to minimal velocity estimates of Sigal and Soffer ||58|| , which were proved using the 
ideas of Mourre [53]; see also Perry, Sigal & Simon | JS|. For our application to the the operator 



e g&(B), we shall refer to special cases of results stated in Skibsted |59| and Debievre, Hislop 



& Sigal 0. 

We first introduce some definitions and assumptions:. 
(AO) 

(a) Let N > 2 and g& G C^°, be a smooth characteristic function which is equal to one on 
the open interval A. 

(b) Let H and A denote self adjoint operators on a Hilbert space 7i. Assume H is bounded 
below. Let T> denote the domain of A and T>(H) the domain of H. Assume T> H T>(H) is dense 
in V(H), and let (A) = (I + |A| 2 )i 

(Al) Let ad^(-ff) = H. For 1 < k < N*, define iteratively the commutator form 



i k ad k A (H) 



„-fc-i„ jft-i 



ad«C (H) , A 



(2.24) 



oeOH T)(H). Assume that for 1 < k < A^*, i fc ad^(if) extends to a symmetric operator with 



domain T>(H), where Af* 



N + 



+ 1. 



(A2) For \s\ < 1, e' lAs : V(H) -> V(H) and sup 



|s|<l 



He' lAs ijj \\ n < oo, for any ip e T>(H). 
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(A3) Mourre estimate: 

g A {H) i[H, A] g A (H) > 9 g A (H) 2 (2.25) 

for some 9 > 0. 

Under the above assumptions, we have, via Theorem 2.4 of [pSH , the following: 
Theorem 2.4. Assume conditions (A0)-(A3). Then, for all S\ > and t > 

|| F f4 < 9A(ff)(^>^ f ^ ll« < C(t)-^i || ^ (2.26) 



t 

where 9 is as in ( 2.2dQ . 

To prove this theorem we set A(r) = A — br, where r = t + 1. Then, for any b < 9, we 
have by ( |2.25| ) that the Heisenberg derivative, £M(t) = <9iv4(r) + A(r)], satisfies 

9aDA(t) 9a = g A (z[H,A}-b)g A > (9 - b) g\. (2.27) 

The result now follows by an application of Theorem 2.4 of pD[ . 

Remark: The hypotheses of Theorem 2.4 can be relaxed |2D|] to the following two conditions: 
(i) the operator g A (H) ad\(H) g A (H) can be extended to a bounded operator on 7i for n = 



0,1,2..., 



2 ~ 



+ 1, and (ii) the Mourre estimate, Q2.25|) . 

We shall apply Theorem 2.4 for the case H = B = ^ — A + m 2 + V(x), A = (x ■ p + p ■ x) /2, 
p = —iV x and 7i = L 2 . Before verifying the hypotheses of Theorem 2.4, we show how this the- 
orem is used to derive the desired estimates on e tBt g A (B). With the application of Theorem 2.4 
in mind, we estimate the norm of (x)~ a e tBt g A (B) by decomposing it into the spectral sets on 
which A/t is less than and greater than or equal to 9. Let Ax denote an interval containing A 
and denote by C the operator that associates to a function / its complex conjugate /. Further- 
more, note that e tBt = C e~ lBt C and C g{B) = g(B) C, if g is real-valued on the spectrum 
of B. Then we find, 

II (x)" 7 e iBt g A {B)^ || 2 = || (x)-° g Al (B) e lBt g A (B)ip || 2 

iBt 



(x)-°g Al (B)C e- lBt C g A (B)^\\ 2 



(x)-°g Al (B) (Ay- (A)- e- iBt g A (B)ip 



2 



< || {x)-°g Al (B) (Ay\\ B{L . } ■ || (A)- e- wt g A (B)ip || 2 

< C || (A)-° e- wt g A {B)i, || 2 

< d|| F (j < 9^j (A)- e- wt g A (B)^p || 2 

+ C 2 \\ ^ (y > # ) (A)'* e- wt g A (B)%p || 2 (2.28) 
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We have used here that, with A defined as above, the L 2 operator norm of {x)~ a g Al (B) (A) a 
is bounded; see [ I8| . By Theorem 2.4, the first term on the right hand side of ( |2.28j ) can be 
estimated from above by: 



.4 



d|| F < 6) {A)- e~ im g A (B)g Al (B)^ || 2 < 
C.Citr^W (A)Tg Ai (B)(x)-%\\ B(L2) \\{x)^\\ 2 . 



Since 9t > 0, the second term is bounded as follows: 



(A)- e- wt g A (B)ip || 2 < C(t)^ 



From (ggp , (|2~29| ) and (^30D , we have 

||(a;}- CT e^ A (S)V||2 < C (<*>-» + (t)'-f + 
Use of ( p.31| ) in a direct estimation of ( |2.19| ) gives: 

Uxyosmwhw^c ((t)-f +2+ - + <tr+ 2 ). 

We make the final choice of N and o after estimation of S\{t\ 



(2.29) 



(2.30) 



(2.31) 



(2.32) 



To complete the estimation of Sl(t), it remains to verify hypotheses (AO)- (A3). These 
hypotheses are known to hold for the operator H = B 2 = — A + V + m 2 ; see ||18|| . Hypothesis 
(AO) holds H = B with domain W /1 ' 2 (lR n ). Clearly hypothesis (A2) holds for H = B since it 
holds for H = B 2 . Hypothesis (Al) can be reduced to its verification for B 2 using the Kato 
square root formula ||49|| : for any ip G T>{B 2 ): 



B^ = 7T- 1 / W - l/2 B 2 (B 2 + w)- 1 t/j dw. 
Jo 



(2.33) 



To verify (Al) for the case k = 1 we must show that [A, B] B~ x is a bounded operator on 
L 2 . Using ( |2.33| ) we have: 

[A, B] B~ l = n- 1 / w* (B 2 + w)- 1 [A, B 2 } B' 1 (B 2 + w)' 1 dw (2.34) 
Jo 

~B 2 ,iA\ = 2B 2 -x-VV -2{V + m 2 ) , (2.35) 

[B 2 ,iA}B- 1 = 2B - (x-VV + 2V) B~ l - 2m 2 B-\ (2.36) 



Since 



we have that 
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Substitution into Q2.34Q we get: 

9 roo 

~ / 

7T JO 



2 r°° i 

[A,B]B~ 1 = - / to* B (B 2 + w)- 2 dw 

IX Jo 

- - / (B 2 + w)- 1 (x- W +2V -2m 2 )B- 1 (B 2 + W)- 1 dw 

it Jo v ' 

= Jx + J 2 . (2.37) 
The term J\ can be rewritten as 

2 t°° i d 

Ji = / — (B 2 + w)- 1 dw B = I, (2.38) 

7T JO dw 



which follows from integration by parts and the formula ( [2.33D with ip replaced by B ip. Using 
that IK^ + w)- 1 !! < (Q 2 + w)-\ 



|| (x-W + 2V - 2m 2 ) B- 1 || < \\x-VV + 2V - 2m 2 \\ 00 WB^Ws^, (2.39) 

and the hypotheses on V we have that for some constant C dependent on V, 

|J 2 | < C w* {Sl 2 + wY 2 dw < oo. (2.40) 
Jo 

The higher order commutators are handled in a similar manner; they are even simpler because 
each successive commutator results in an extra factor of (B 2 + w) . 

This leaves us with verification of the Mourre estimate (A3) for H = B. Proposition 2.2 of 
20H says that under hypotheses on B 2 , verified in |18j, that (A3) holds for B as well. 



Estimation of S2 : 

In 5*2 (t), the energy is localized away from A so we seek to use the W k,s estimates of Propo- 
sition 2.1. Note that of the cases I = 1 and / = 2, the case / = 1 is "worse" because g A localizes 
the energy away from the singularity at A, and the I = 1 term has slower decay for large energy. 
We therefore carry out the estimation for 1 = 1. 

Let q~~ l + p'^ 1 = 1/2 and p'^ 1 + p^ 1 = 1. Since a > |, (x)~ CT G L q and therefore: 

||(x)~ CT Sm (x}-° f\\ 2 < C\\(x}-°\\ B{L ^ L2) \\e Wot (B -A + zO)- 1 g A (B ) W* + (x)-° f\\ p , 

< C\\e iB ^B^\\ B{w ^ Lpl) Wg^Bo^Bo-A + tOy'BoWKx)-* f\\ 1>p 

< CWe^B, 1 ^^^ \\g^B Q ){B Q -k + iQ)- l B,-B Q Wl{x)-° f\\ p 

< CWe^B, 1 ^^^^ \\g A (B )(B -A + iO)- l B \\ Bm \\B W;(x)-° f\\ p (2.41) 



The three factors in ( |2.41| ) are estimated as follows: 
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(i) By part (a) of Theorem 2.1, 

W^B^hiw^Lp') < C\t\-^2, (2.42) 



where p = 2 - 8/(n + 6) and p' = 2(n + 2)/(n - 2). 

(ii) ||<7 A (-B )(-B — A + iO)^ 1 B \\ B ( LP ) is bounded because g A (/x)(/z — A + zO) _1 /i is a multiplier 
on L p 0. 

Using the boundedness of on W 2 ' p and that a > max{§, 2}, we have 

II b wi(x}-° f\\ p < \\w;(x}-° f\\ ltP 

< \\W*\\ B(W 2, P) ||(x)- ff /|| 1>p 



< C\\J II 

Therefore, 



||(x)-^|(x)->|| 2 < |^||i, 2 . (2.43) 

We now combine the estimates of Sj(t), j — 1,2 to complete the proof. Combining ( p.32j ) 
and ( |2.43| ) and taking e j yields 

IKa;)-' 7 e iBt (5 - A + i0)~ l P c (x)- a tjj\\ 2 < C ((i)-f+ 2 +^ + ( t }-+ 2 + (t)-£s ) 



.1,2- 

(2.44) 

The estimates ( |2.17| ) now follows by taking A" and a such that 

2n 

a > 2 + 

n + 2 

n 

a > max{— ,2} 

N 2n 

— > 2 + + Ei 

2 ~ n + 2 

The constraints on o are implied by the hypothesis a > a*(n). The remaining constraint holds 
if AT > 8. Since we have applied Theorem 2.4 in our estimation of S[ , we need that (Al) hold 



with N* 



1 > 10. 



This completes the proof of Proposition 2.2. 



3. Existence theory 

In this section we outline an existence theory for (1.1) with initial conditions u(x, 0) = Uq G 
Vy 2 ' 2 (lR 3 ) and d t u(x, 0) = u\ G W /1 ' 2 (1R 3 ). We first reformulate the initial value problem and 
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then introduce the hypotheses on the operator H. Regarding ( |1.1|) as a perturbation of a linear 
constant coefficient equation, we first write the initial value problem as: 

-Vu + Xf(u). (3.1) 

u (x) 

cos(B t) sm(B t)/B \ 
B sm(B t) cos(B t) J ' 

-Vu + Xf(u)J' 

Then, the initial value problem can be reformulated as a system of first order equations: 

d t u = (_^ 2 J)u + F(u). (3.2) 

We now follow the strategy of reformulating the problem of finding a solution of the initial 
value problem as the problem of finding a fixed point u of an appropriate mapping. In particular, 
we seek u in the space 

X = ^(IR 3 ) x W^' 2 (1R 3 ) (3.3) 

satisfying 

iu = u, (3.4) 

where 

Au(t) = U(t)u + [ U(t- s) F(u) ds. (3.5) 

Jo 



d 2 t u + Blu = 

u(x, 0) = 

d t u(x, 0) = 

Let 




For a discussion of the existence and low energy scattering for the case V = 0, see [p3[ , |p5| , p6 | 
and references cited therein. 

Hypotheses of H = —A + V 

(Hi) v e W 1 ' 00 

(H2) H = B 2 is a positive and self-adjoint operator 

(H3) The semi-infinite interval, [m 2 , oo), consists of absolutely continuous spectrum of H . 
(H4) H has exactly one (simple) eigenvalue Q 2 satisfying < Q 2 < m 2 , with corresponding 
eigenfunction ip e L 2 , \\fW2 = 1- 
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(N) f(u) = u 3 + fi(u), fiiu) = 0(u 4 ) and f(u) is smooth in a neighborhood of u = 0. 

These hypotheses are by no means the least stringent, but are sufficient for the present 
purposes. 

Theorem 3.1. ("Local existence theory,) 



Consider the Cauchy problem for the nonlinear Klein Gordon equation (LI) with initial data, 
u of class X . 

(a) There exists strictly positive numbers, T max and T max (forward and backward maximal 
times of existence) which depend on the Xo norm of the initial data, such that the initial value 
problem has a unique solution of class C° ((— T max , +T max ); X ) in the sense of the integral 
equation (\3.4 )- 

(b) Fort G (— T max , +T max ) conservation of energy holds: £[u(- , t) , d t u(- , t)] = £[uq,ui], where 
£ is defined in equation fll-4 )- 

(c) Either T max is finite or T max is infinite. If T max is finite, then 

l| m H-,t) ||xo = oo. (3.6) 

* | L max 

The analogous statement holds with T max replaced by T max and linu^T^ replaced by lim^-^™^ 



We now sketch a proof of Theorem 3.1. We restrict our attention to the case t > 0; the 
proof is identical for t < 0. Furthermore, we shall, for simplicity, consider the case of the cubic 
nonlinearity: f(u) = u 3 . Using hypothesis (HI) on V, it is simple to show, using the estimates 
of Corollary 2.1, that for T > sufficiently small and depending essentially on the Xo norm of 
uo, the mapping A maps a closed ball in C° ([0, T); Xo) to itself, and is a strict contraction. This 
ensures the existence of a unique fixed point, which is a C°([0, T);X ) solution of the initial 
value problem ( [3.1|) in the sense of the integral equation (|3.4j) . That the Xo norm must blow up 
if T max is finite is a standard continuation argument based on the fixed point proof of part (a). 

In section 7, we study the asymptotic behavior of solutions as t — > ±oo. A consequence of 
this analysis and the argument below is that if u satisfies the more stringent hypothesis that 
the norm ||uo||x is sufficiently small, where 

x = (w 2 ' 2 n w 2 ' 1 ^ x (w 1 ' 2 nw 1 ' 1 ) , (3.7) 

then T max = T max = oo and the solution decays to zero as dilineated in the statement of 
Theorem 1.1. 
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We reason as follows. A unique solution, u(t), exists locally in time and is continuous in t with 
values in Xo- It follows that for \t\ < T, the W 1 ' and L 8 norms of u(-, t) are finite. Moreover, the 
solution satisfies energy conservation and therefore ||p^i.2 + ||f(t)||2 is bounded uniformly 
by a constant which is independent of T, and depends only on HwoIIvp 1 - 2 + Ikolh- To ensure 
that the quantities estimated continue to be well defined and satisfy the 'a priori estimates of 
section 7 it suffices to show that ||u(t)||x remains bounded as t \ T. A direct and elementary 
estimation of the integral equation fl3.4|) yields the estimate: 



|u(t)||x < C(lko||w 2 . 2 + IKIIw 1 . 3 ) 

+ C{\\V\\ w uoo,\) f (||«(s)|| w i, a + \\u(s)\\^ a + \\u 2 (s)Vu(s)\\ 2 )ds. 

(3i 



The main tool used in obtaining (|3.8|) is the bound \\E[ \t)di\\j3(L 2 ) < 



m 



A simple consequence of the 'a priori bound (|7.52|) and the decomposition of u(t, ■) is that 



\\Vu{t)\\ 4 + \\u{t)\\ s < C, \t\<T. (3.9) 



The constant C, depends on the norm ||u ||x- The first two terms in the integrand of ( |3.8| ) 
are uniformly bounded by conservation of energy. Furthermore, this energy bound together 
with ( |3.9| ) implies a bound on the last term in the integrand of the estimate ( |3.8| ). It follows 
that ||u(t)||x remains bounded as t f T. Therefore, given the estimates of section 7 we have 
Tmax — °° > an d the decay of solutions. 
A further consequence of the proof is: 

Corollary 3.1. Under the hypotheses of Theorem 1.1, the solution of the initial value problem 
exists globally in W 2 ' 2 (JEl 3 ) and satisfies the estimate: 

\\u(t)\\ W 2,2 < C (t). (3.10) 



4. Isolation of the key resonances and formulation as coupled finite 
and infinite dimensional dynamical system 

Using the notation ( |1.5[ ), the initial value problem for ( |1 . 1| ) can be rewritten as 

d 2 t u + B 2 u = A f(u), f(u) = u 3 (4.1) 
u(x,0) = uq(x), dt u(x,0) = ui(x). 
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For small amplitude solutions, it is natural to decompose the solution as follows: 



u(x,t) = a(t)(p(x) +rj(x,t), (4.2) 
{r){-,t),(p) = for all t . (4.3) 



Substitution of ( gig ) into ( P~ip gives 

a"(p + d 2 r] + n 2 a<p + B 2 r] = X f(a(p + rf) (4.4) 
We now implement ( |4.3|) . Taking the inner product of Q4.4|) with <£> gives 

a" + fi 2 a = X((p,f(a(p + ri)). (4.5) 
Let P c denote the projection onto the continuous spectral part of B 2 , i.e. 

P c v = v — (</?, i>)<£>. (4.6) 

Then, since ?7 = P c ?7, we have 

+ B 2 V = \PJ(a<p + v) (4-7) 



Equations ( |4.5| )- (|4~7| ) comprise a coupled dynamical system for the bound state and contin- 
uous spectral components (relative to H = B 2 ) of the solution u. The initial conditions for this 
system are given by ( |1.16| ). 

Expansion of the cubic terms in ( f4.5| - [4.7D gives the system 



a" + Q 2 a = X 



(4- 



a 3 J V + 3a 2 j ^r] + 3aJ (p 2 rj 2 + / W 
d\ i] + B 2 i] = AP C (afy 3 + 3a 2 (p 2 ij + Zatpr] 2 + r/ 3 ) , (4.9) 

with initial conditions 

o(0) = (tp,Uo) , o'(0) = 

7/(35,0) = P C M , 3t77(M) = Pc«l (4-10) 



We now locate the source of the resonance. Equation ( |4.8| ) has a homogeneous solution 
which oscillates with frequencies ±f2. Thus, to leading order, 77 solves a driven wave equation 
containing the driving frequencies ±3f2. If 9Q 2 > m 2 , then by (H3), we expect a resonant 
interaction with radiation modes of energy 3f2 e a con t(H). 
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This resonant part of T], has its dominant effect on the a-oscillator in the term of (|4.8|) , which 
is linear in 77. Our goal is to derive from (|4.8| - |4.9|) an equivalent dynamical system which is of 
the type described in the introduction, but which is corrected by terms which decay sufficiently 
rapidly with time and which can therefore be treated perturbatively. 

We first write 

V = m + V2 + V3 , (4.11) 
where i]i(t) satisfies the linear dynamics with the given initial data: 

d 2 Vl + B\ = 0, m(x,a) = P c u , d t V i(x,0) = P cUl , (4.12) 

and f]2{t) is the leading order response 

8 2 t 772 + B% = X a 3 P c <^ 3 ; r] 2 {x, 0) = 0, d t r] 2 (x, 0) = . (4.13) 

Equation (^4.8| ) can now be written in an expanded form. 



a" + tt 2 a 



X 



J V A + 3a 2 j ^ 3 ( + 772 + 77 3 ) + 3a J V 2 n 2 + J W 3 



= F{a,rj) 



(4.14) 



We expect the function a(t) to consist of "fast oscillations", coming from the natural frequency 
Q and its nonlinearly generated harmonics, and slow variations due to its small amplitude. We 
next extract from a(t) the dominant "fast" oscillations of frequency Q: 

A e int + A e~ m . 



a(t) 



We then substitute (|4.15|) into ( |4.14|) and impose the constraint: 

A , e iat + A'e- int = 
Equation (4.14) then is reduced to the first order equation: 



A' = (2zO) 



F(a, V ) 



(4.15) 



(4.16) 



(4.17) 



where F(a,r)) = F(A, A,r),t). From ( |4.14j ) we have that F(a,r)) is the sum of the following 
terms: 

Fi(a) = 
F 2 ( y a,r ]j ) -- 
F 3 (a,r]) = 

FM = X [ iprf . (4.18) 



Xa" J ^ 

3Xa 2 J <p\ , J = 1,2,3 
3Xa J (f 2 rj 2 
A / W ■ 
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The remainder of this section is primarily devoted to an (involved) expansion of F 2 (a,r) 2 ). 
The terms of this expansion are of several types: (a) the resonant damping term, (b) terms of 
the same order whose net effect is a nonlinear phase correction and (c) higher order terms which 
are to be treated perturbatively in the asymptotic analysis as t — > ±00. 

Computation of F 2 (a, 772) 
We first break 772 into a part containing the key resonance rf 2 and nonresonant part r\ 2 r . 



Proposition 4.1. 



where 



and 



+ 



02 



V2 = rf 2 + rfi 



—e lBt C e-^ B - 3 ^A 3 (s) ds P c <p 3 
2iB Jo \ > w 



JBt 



2iB 



3 /V is(B - n U 2 (sU(s)cfe 
. Jo 

+ 3 [' e- is{B+n) A\s)A(s) ds 
Jo 

+ [ t e-^ B+3 ^A\s)ds] P c <p 3 
Jo J 

/■* e «(B+3n) A 3 (s) ds 
Jo 



-iBt 



2iB 



+ 3 /* e is{B+n) A 2 (s)A(s) ds + 3 C e ls{B ~ Q) A\s)A(s) ds 
Jo Jo 



[ e is( - B ~ 3 ^ A 3 (s) ds 
Jo 



Pc^ 3 



3=1 



(4.19) 



(4.20) 



The superscripts r and nr denote respectively a resonant contribution and nonresonant contri- 
bution. 



Proof: The solution of ( |4.13j ) can be expressed, using the variation of constants formula, as: 

sin B(t — s 



r)2 = A / V D — -a 3 (s) ds Pcv? 3 , 

70 B 



(4.21) 



Substitution of ( |4.15| ) for a(s) and using the expansion of sin (B(t — s)) in terms of the operators 
exp (±iB(t — s)) leads to an expression for rj 2 which is a sum of eight terms. The term 7/2, as 
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defined above, is anticipated to be the most important. The other seven terms are lumped 
together in 7/2 r . 

We now focus on r) 2 . In order to study rf 2 near the resonant point 3Q in the continuous 
spectrum of B, we first introduce a regularization of T] 2 . For e > 0, let 



V 

V2s 



A Wt f e -is(B-3Q+ie) A 3^ dg p^3 



, yp , . r (4-22) 

2iB Jo 

Note that rf 2 = lim e ^ ?72 £ . The following result, proved using integration by parts, isolates the 
key (local in t) resonant term. 



Proposition 4.2. For e > 0, 



V 

V2e 



^ [B(B - 3fi + ie)] 1 e 3ff *A\t) e £t P c cp 3 
A 



B(B — 3fi + ie) 



_1 JBt 



Pc^ 3 



- -A 
2 



B(B - 3tt + ie)] e iBt / e -<«(B-3n+fe)^/ rfs p^s 



vis + C 1 + < 2 



(4.23) 



Remark: The choice +ie in (|4.22|) is motivated by the fact that the operator 



(4.24) 



satisfies appropriate decay estimates as t — > 00; see Proposition 2.1. It follows that the limit as 
e — > + of the last two terms in ( |4.23|) decay in time like (t) _1_a (« > 0); see section 7. 

We now use the above computation to obtain an expression for F 2 (a, rj 2 ). First, from Propo- 
sition 4.1 we have 



F 2 (a,r] 2 ) 



F 2 (a, V r 2 ) + F 2 (a,^ r ) 

limF 2 (a,C) + limF 2 (a,C 1 + C 2 ) + F 2 (a,^) 



F 2 (a, V :) + F 2 (a,C rl + vT) + F 2 (a }V n 2 



(4.25) 



where F 2 (a, •) is defined in (|4.18 ). We begin with the contribution to ( 4.17[) coming from 
F 2 (a,i]l). What follows now is a detailed expansion of the term F 2 (a,r]l) and F 2 (a,r) 2 r ). The 
terms F 2 (a, rf^ rl + rf^ r2 ) can be treated perturbatively by estimation of its magnitude; see section 
5. 

Computation of F 2 (a,rjl) 
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Let 



A 



.. /_ 3 1 B-3Q 



(4.26) 



r = 



71 
7T 



3Q 



Pc ^ 



(P c y> 3 ,5(B-3f2)P c ^ 



(4.27) 



By hypothesis Q, T > 0. 

We now substitute the expression for r]l £ , given in ( |4.23[ ) into the definition of F 2 (a, r)l e ) in 
( |4.18D . Passage to the limit, e — > 0, and use of the distributional identity: 



(x ± zO) 



1 = lim (x ± ie) 1 



P.V. x =F iirS(x) 



(4.28) 



yields: 



Proposition 4.3. 



limF 2 (a,7/ 



-a 2 (a-*t) 



A\ 4 Ae im + A b e blUt + 2\A\'A 6 e 



5„5ifit 



1 2 A3„3iQ,t 



(4.29) 



We have completed the evaluation of the first term in Q4.25| ). To calculate the contribution 
of ( |4.29| ) to the amplitude equation ( |4.17| ) we need only multiply (|4.29|) by (2zf2) _1 e~*™. This 



gives 



Proposition 4.4. 

(21Q)- 1 e- im F 2 (a, V : 



3A^ 



(iA + r) 



A\*A + A 5 e 4int + 2\A[ z A*e 



2 A3„2iQt 



(4.30) 



The term — f tt- r|A| 4 A plays the role of a nonlinear damping; it drives the decay of A and, in 



turn, that of rj; see the discussion in the introduction. 

Computation of F 2 (a, 7/2 r ) • 
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We now focus on i^a, rf^), the third term in ( |4.25|) . This requires a rather extensive, 
expansion of rf^ = J2jV2j'i see ( ^-20|) . From Proposition 4.3, we expect the dominant terms 
to be (9(|A| 5 ). Our approach is now to make explicit all terms which are formally (9(|t4| 5 ) 
(anticipating that A' = 0(\A\ 3 ), A" = 0{\A\ 5 ) and {\A\ 2 )' = 0{\A\ 6 ) ) and to treat the 
remainder as a perturbation which we shall later estimate to be of higher order. This expansion 
of rff is presented in the following proposition which we prove using repeated integration by 
parts. As written, these expressions are formal. By the definition of i^a,?^) we require that 
they hold when integrating against a rapidly decaying function, i.e. if 3 . 



Proposition 4.5. The following expansions hold in S': 

A , . n . „+r> - q 3 A 



+ 



2B(B-tt) 

3Xe iBt 
25(5-0) 
3A 

2iB(B - O) 2 



\A\ 2 Ae itn P c <p 3 + 



2iB(B - O) 2 
A \ 2 A + (\A\ 2 A)'\ t=0 ] P c p 3 

e iBt f* e -is{B-n) n A \2 A y ds p 
Jo 



e im {\A\ 2 A)' P c <p 3 



(4.31) 



3A \A\ 2 ~Ar- im V in 3 I 3A 

2B(B + n) 11 c ^ + 2iB(B + o) 2 



e- ltn (A 2 A' + 2A'\A\ 2 )P c ip 3 



+ e 



iBt 



A 



\A n \ 2 A- 



3A 



2B(B + 0) 



2 2B{B + 0) 

3A e* Bt /* e- ts{B+n \\A\ 3 A)"(s) ds P cV 
Jo 



(\A\ 2 A)'\ t=0 P cV 



2iB(B + O) 2 



(4.32) 



V23 



\A(tf 



sun 



2B(B + 30) 
3^ e iBt 

2 5(5 + 30) 

3 A 

2 i5(5 + 30) : 



Pc¥ 3 + 



3A 



2i5(5 + 30) 



r e- 4 ' 3n (A(t) 3 )' P c <^ 3 



^ + (A 3 )'\ t=0 P cV 3 

(A(s) 2 )" ds P cV 3 



_ e iBt j " e -is(B+3n) 



(4.33) 



A 



25(5 + 30) 



3A 



e 3itn (A 3 )' P c r 



2iB(B + 30) : 
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-iBt 



+ 



2 5(5 + 30) 
3A 

2i5(5 + 30) 



^0 + (^ 3 )'l t =o]Pc^ 
-iBt /V< B4 ^>(A 3 (s))''ds p< 



(4.34) 



3A 



\A\ 2 A e im P c ^ a - 



25(5 + 0) 

o -iBt 
-A 



3A 



2^5(5 + o) 2 



e «n(|^|2^)/ Pc¥ ,3 



+ 



2 5(5 + 0) 
3A 

2i5(5 + O) 



\A \ 2 A + (\A\ 2 A)'\ t=0 \P c ^ 

-iBt f* e is(B + n)n A \2 A y, dg p a 
2 Jo 



(4.35) 



^26 



3A 



25(5-0) 

o -iBt 
-A 



\A\ 2 A e~ im P c <p* - 



3A 



2iB(B - O) 



-e~ im {\A\ 2 A)' P c <^ 3 



+ lx 



2 B(B-fl) 
1 



|A |% + (|A|^4)'| t=0 P c ^ 



2i 5(5 -O) 2 



e iBt / e u(B-a), 



(\A\ 2 A)" ds P c ^ 



(4.36) 



AA 3 (t) e 



%7 



25(5- 30 -iO) 
3^ e~ iBt 
2 5(5 - 30 - iO) 
3A 



3A 



2i5(5-30-iO) : 



+ 2,5(5- 30 -iOr '^Jo ^ 
Recall that our goal is to elucidate the structure of the amplitude equation: A' = {^lifl)^ 1 e~ int F , 
in ( [4.171 ), where we first focused on the contribution: (2ifl)~ 1 e~ int F 2 . From Q4.25Q and Propo- 
sition 4.4 we see now that we need to obtain convenient expressions for 



(2tfl)- 1 e~ iQt F 2 (a,r ] ™) 



^(2iO)- 1 e- int F 2 (a,^J) 

3=1 

7 r 

3\{2tfl)- l a 2 e- tnt J2 I vSZj ■ 



(4.38) 



3=1' 
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Each of the seven terms contributing to rf^ is expressed as a part which is (9(|t4| 5 ) plus an 
error term which is estimated in magnitude in section 5. The (9(|t4| 5 ) part and error terms are 
displayed in the following two propositions. 



Proposition 4.6. Let 

Then, 



p(C) = (Pc^,^ 1 (5-C)- 1 Pc^ 3 )- 



(4.39) 



(2ifi)- 1 e- iS2t F 2 (a,r ?2 "[) 

£ p (Q) \h A \ 2 A 3 e 2int + -\A\ 4 A + hA\ 4 Ae- 2iQt 
a2 4$ 2^ 4z 



+ E nr 



21 



(4.40) 



(2iQ)- 1 e- iQt F 2 (a,r ] ^) 
= ^r-p(-fi) \^\A\ 4 A + ^\A\ 4 Ae- 2im + ^\A\ 2 A 3 e- 4tnt ] + E% (4.41) 



(^-^-^(a,^) 

= ^-p(-3fi) [i A\ 4 Ae~ 2int + -\A\ 2 A 3 e- 4tnt + ^V 6 ^'] + £ 2 n 3 r (4.42) 
12 ~ 4?- 2 4 



(2ifi) 
A 2 



-1 -ifit 



^2 MI) 



p(-3fi) 



+ ^ 3 |A| 2 e 2in * + ||A| 4 A] + EH 



(4.43) 



(2in)- 1 e- in *F 2 (a,i^) 

= ^r-p(-fi) \-\A\ 2 A 3 e 2tm + -\A\ 4 A + -\A\ 4 Ae- 2int ] + E™ (4.44) 
12 4^ 2^ 4z J 



(2ifi) _ 
A 2 



-int 



r9 



ft' 



9 



p(Q.) ^\A\ 4 A+^\A\ 4 Ae 



4"X -2ifit 



4i 



2i' 



+ i|A| 2 A 3 e- 41 ™ 
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(4.45) 
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(2tn)- 1 e- iQt F 2 (a,ri^] 



A 



3 , 



3 , 



— p(3ft + iO) -\A\ A Ae~ 2im + -\A\ 2 A 3 e- AiQt + A 5 e~ 6int 



rpnr 



(4.46) 



In our analysis of the large time behavior (t — > ±oo), we shall require an upper bound of 
the error terms EVfj given in: 

Proposition 4.7. 

\E% | < C„ |A| 2 \A\ 2 { \A\ 2 \K\ + IK^^S-Or 1 e ifl * P c <p 3 || 2 
+ \\(x)-*(B-Cj)- 2 fe iB <t-*) 0((\A\ 3 )") dsP c cp% }, (4.47) 

« 

where Ci = C6 = ^ C2 = Cs = ~~ ^> C3 = C4 = — 3£7, and ( 7 = 3fl + iO. 

In section 7, we shall estimate this expression using the decay estimates of section 2, in particular 
Proposition 2.2. 

The desired form of the A-equation is now emerging. Use of Propositions 4.4 and 4.6 in 
(p7|) yields: 

Proposition 4.8. The amplitude A(t) satisfies the equation (see Proposition 4.6 for the defi- 
nition of p(C))'- 



A' 



-iX 



2tt 



I i 1 ■:> ||2 1 \Z.XtiH ■:>: ||2~T . _|_ ^ 3 g-4int 



3\A\ 2 A + A 3 e 2im + 3\A\ 2 Ae~ 



— r\A\ 4 A - — \ 2 
4 n 11 4fi A 

3A 2 



A - 5p(fi) + 3p(-Q) + p(-3fi) 



5 4ifM 



A° e 



2^3 2«fit 



\A\ Z A 6 e 



4ft 
3?A 2 

' 4ft 
3? A 2 

■ — — I A I 4 "J <=- 2i ™ 

4ft {A{Ae 

3iA 2 , .,0^ 



zA + T - p(-3ft) 



4ft 



A ,2 Z 3 e _ 4im 



- 2 A + 2ir + p(ft) + ip(-3ft) + 3p(-ft) 
7p(ft) + 9p(-ft) + p(-3ft) + p(3ft + i0) 
3p(-ft) - 2zp(-3ft) + p(3ft + i0) + 3p(ft) 



3iA 2 _ 5 



4ft 



A e 



ip(-3ft) ip(3ft + z0) + E 

3 



(4.48 



32 



where 



E = (2in)- 1 e int 3Xa J y?r/ 2 
+ (2iQ)- 1 \e- int 



nr 
■2.i 



+ (2 i Q)- 1 e- int [F 2 (a,r ]l ) + F 2 (a, V3 ) 
+ F 2 {aX rl + € r2 ) 



(4.49) 



Here, F 2 (a, rjj) is given by ( |4. . 

In the next section we show how to rewrite ( 4.48| ) in a manner which makes explicit which 
terms determine the large time behavior of the amplitude and phase of A(t). 



5. Dispersive Hamiltonian normal form 

To analyze the asymptotic behavior of A(t) (or equivalently a(t)) and f]{t,x) as t — > oo it is 
useful to use the idea of normal forms ||, p6| , |52| from dynamical systems theory. We derive a 
perturbed normal form which makes the anticipated large time behavior of solutions transparent. 



Proposition 5.1. There exists a smooth near-identity change of variables, A i— ► A with the 
following properties: 

A = A + h(A,t) 
h(A,t) = 0{\A\ 3 ), \A\^0 

h(A,t) = h(A,t + 2nQ- 1 ), (5.1) 



and such that in terms of A equation ( 4.48 ) becomes: 

A' = i\c 21 \A\ 2 A + X 2 d- i2 \A\ 4 A + i\ 2 c 2>2 \A\ A A 

+ 0(|i| 7 ) + E, (5.2) 

where X 2 d^ 2 = — §^~r < 0. The constants c 2 \ and c^ 2 are real numbers, explicitly calculable in 
terms of the coefficients appearing in ( |4.4<Sp . The remainder term |E| is estimable in terms of 
|E| for \ A\ < 1 (equivalently \A\ < 1). 
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Remarks: 

The point of this proposition is that in the new variables the dynamics evidently have a dissi- 
pative aspect. Specifically, neglecting the perturbation to the normal form one has: 



a \a\ 



3 A 2 , ~. 

-— — — r \a\ 

8 n 11 



< o 



(5.3) 



We therefore refer to ( |5.2|) as a dispersive Hamiltonian normal form. In finite dimensional 
Hamiltonian systems, the normal form associated with a one degree of freedom Hamiltonian 
system is an equation like (|5.2|) , but with all the coefficients of the terms |/1| 2 "M being purely 
imaginary. Here we find that resonant coupling to an infinite dimensional dispersive wave field 
can lead to a normal form with general complex coefficients, which in our context implies the 
internal damping effect described above. See section 8 for further discussion. 

We now present an elementary derivation of the change of variables ( |5.1|) leading to (|5.2|) . 
Equation ( }4.48|) is of the form: 

A> = E E a ki Ak 3* e i{k - l - l)m + E, (5.4) 

j€{3,5} k+l=j 



where the coefficients a^i can be read off ( f4.48|) . The proof we present is quite general and 
shows, in particular, that the normal form for equations like (|5.4|) is (|5.2|) 



Proposition 5.2. There is a change of variables, as in ( l5.Jp , such that equation fl5-4 ) is mapped 
to: 



A' 



= k 2l \A\ 2 A + k m \A\ A A 
+ 0{\A\ 7 ) + E , where 



(5.5) 



E = Eo(/ + /i), 
hi = «2i 

h2 = « 32 + (2iQ)- 1 



2 |C * 031 



2a 30 a 12 + 2\a 12 \ 



(5.6) 



The conclusion in Proposition 5.1 concerning the "damping coefficient" c/32 depends on the 
particular properties of the coefficients in (|4.17|) . In particular we have from (|4.17| ) and (|5.4| ) 
that: 

3A 



«21 


= a 12 = 




A 


«03 


2iQ^ 






«32 


4Q 



2iQ 
i 



I Il4 



MX 2 

An 



[A - Bp(Q) + 3p(-0) + p(-3fl) 



(5.7) 
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From these formulae, we have X 2 d 32 is given by the real part of a 32 , the "singular" Fermi golden 
rule contribution. 

Remark: The construction of the map A i— > A can be applied as well to equations of the form 
( |5.4| ) where the right hand side is an arbitrary expansion in powers of A and A. 

To make the structure clear we write out the equation with a particular ordering of terms: 

A' = a 2 i \A\ 2 A + a 32 \A\ A A + 3 {A) + 5 {A) + E, (5.8) 
where 

3 (A) = a 30 A 3 e 2int + a 12 AA 2 e' 2int + a 03 Z 3 e~ 4int , and (5.9) 
O s (A) = a 5Q A 5 e 4lQt + a 41 A 4 A e 2tm + a 23 A 2 A 3 e~ 2im 

+ a u AA 4 e" 4i ™ + a 05 A 5 e" 6i ™ (5.10) 

Note that each term in 3 (A) and 0§(A) is of the form: oscillatory function of t times order 
0(\A\ 3 ) or 0(\A\ 5 ). 

The computations that follow, though elementary, are rather lengthy so we first outline the 
strategy of our proof. Integration of (|5~8|) gives 

A = A + /* a 21 \A\ 2 A + a 32 \A\ 4 A + 3 {A) + O s {A) + E ds. (5.11) 
Jo 

The idea is that terms with explicit periodic oscillations average to zero and can be neglected in 
determining the large time behavior of the solution. Our strategy is now to expand the explicitly 
oscillatory terms using repeated integrations by parts and to make explicit all terms up to and 
including order \A\ 5 . The computation has two stages. In stage one, after repeated integration 
by parts the equation (|5.11| ) is expressed in the equivalent form: 

A(t) - H!(A(t),t) = A - #i(A ,0) 

+ / "resonant" terms like \A\ 2 A , \A\ 4 A ds 
Jo 

+ / terms of type O^(A) ds + higher order corrections (5.12) 
Jo 

This suggests the change of variables A \— > A\ = A(t) — Hi(A(t),t), giving 

AAt) = A w + [ terms like Ui| 2 v4i , \A 1 \ 4 A 1 ds 
Jo 

+ terms of type 0^{A{) + higher order corrections. (5.13) 

The latter equation is equivalent to a differential equation of the form: 

A' x {t) = terms like \A 1 \ 2 A 1 , |Al| 4 Ai 

+ terms of type O^(Ai) + higher order corrections. (5-14) 
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which is a step closer to the form of the equation we seek. A second iteration of this process 
yields the result. 

We now embark on the details of the proof. 

Expansion of Jq 0^{A) ds: 

Integration of the expression in ( |5.9| ) gives 

-t 



3 {A) ds 



h 3 (A(t),t) - h 3 (A ,0) 



30 / e 2Ms A 2 A , ds + 



+ 



3a- /■' 
~2iVt Jo 

"03 



«12 

2iVt 



f e-^ s (A A 2 )' ds 



e~ ms (A 3 )' ds, 



(5.15) 



where 



h 3 {A,t) 



«30 ^3 e 2iUt _ ^12_ ^ -fit -2ifit _ ^03 -^3 _4ifit 



2iQ 2iQ Aifl 

We now replace A' in ( |5.15| ) by its abbreviated expression given in fl5.8|) . Thus we have 



(5.16) 



3 {A) ds 



h 3 (A(t),t) - h 3 (A ,0) 



+ 
+ 
+ 



2iVt 



«12 

2iVt 

«12 



a 21 \A\ 2 A + 3 (A) + 0(\A\ 5 ) + E 
a 21 \A\ 2 A + 3 (A) + 0(\A\ 5 ) + E 



ds 
ds 



2itt jo 
3a 03 "' 



e -2iUs 2 | A | 

e- 4lQs A 2 



a 21 \A\ 2 A + 3 (A) + 0(\A\ 5 ) + E 



ds 



io . . a 21 |A| 2 A + 3 {A) + 0(\A\ 5 ) + E ds (5.17) 

Substitution of ( |5.9D into (|5.17|) and integrating by parts, we arrive at the following expres- 



sion: 



3 (A)ds = h 3 (A(t),t) + h 3a (A(t),t) -h 3 (A ,0) - h 3a (A ,0) 
o 

1/3 \ /"* 

+ (2 ' a ° 3 ' 2 ~~ 2a3 ° ai2 + 2 l a i2l 2 ) y o l^l 4 ^ 4 rfs 
+ J* 0(\A\ 2 (\A\ 5 + |E| 



(5.18) 



where 



MA*) 



ai2«03 3a 30 a 2 i 



4fi 2 



2«r2t 



3™ 2 

30 /1 5 4iQi 
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+ ( a 2i«i2 + -ao3"i2 - 3a 30 a 3 + 2a 12 a 2 ij A 2 A 3 e 2int 

1/3 \ 
+ ( a i2 + 2 a ° 3 " 21 + 2ai2 "30 J e~ Aim 

+ ^ ( ^03^12 + i« 03 a 3 o ) ^ e~ 6im (5.19) 

Referring back to ( |5.8| ), we see we must now obtain an 

Expansion of J * O^(A) ds: 

From (|5.10|) we have, after integration by parts: 

* 5 (A)ds = h 3b (A(t),t) - h 3b (A ,0) 

+ J* [ \A\\\A\ Z + |E|) ) ds, (5.20) 



where 



hoJA t) — A 5 p 4iQt + — A 4 ~A p 2m - — A 2 ~A 3 p- 2iQt 

_ ^" A A 4 e~ 4 ^ - ^-A 5 e" 6 ^. (5.21) 



Therefore from ( |5.11| ) and our computations we have: 



A - h 3 {A,t) - h 3a {A,t) - h 3b {A,t) 

= A - h 3 (A ,0) - h 3a {A ,0) - h 3b {A ,0) 

+ [ «2i \A\ 2 A + a 32 \A\ 4 A ds 
J o 

+ ( |l«os| 2 - 2a 30 «i2 + 2|a 12 | 2 ) J* \A\ 4 A ds + jf* E ds 

+ J* 0( \A\ 4 (\A\ 3 +|E|)) + o(\A\ 2 (\A\ 5 +|E|)) ds. (5.22) 

This suggests the change of variables: 

Ai = A - Hi(A,t), where 
H x {A,t) = h 3 (A,t) + h 3a (A,t) + h 3b (A,t), (5.23) 

which for small \A\, is a near-identity change of variables. Using this change of variables we 
have that (|5.22|) becomes 

ft 



A 1 = A w + I a 2 i \A l \ 2 A l + a 32 |-4i| 4 ^i ds 
Jo 
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1 

+ 



2iVL 

+ 



3 \ /■* 

-|a 3| 2 - 2a 30 «i2 + 2|«i 2 | 2 J y |Ai| 4 Ai ds 



ill i \ z /JO 

f (2a 2 i |Ai| 2 fTi(Ai,s) +Af Fi(Ai,s) ) rfs + jf* E lC ?s 

+ jT* O ( lAxfflAxI 3 + lExl) ) + O ( lA^d^l 5 + |Ei|)) rfs. (5.24) 
expand the terms involving Hi(A,t). Using integration by parts, as above, we obtain: 
2a 21 jT lAxl^xds =M^i>*) - ha(A 10 ,0) + jf* ( |^-i| 4 (j^-ij 3 + |E X |) ) ds, (5.25) 

jf* A 2 i7 l( fe =/ i5b (A 1 ,t) - /i56(Ai ,0) + jf* 0(|A| 4 (|A| 3 + |E X |) ) rfs, (5.26) 



and 
where 



^21^30 4ifit 4 5 

H 5 (A,t) = h 5a (A,t) + h 5b {A,t) (5.28) 

Thus 



A 1 - if 5 (^i (*),*) = A 10 - # 5 (Ao,0) 

+ / a 21 |Ai| 2 Ai + a 3 2 ote 

JO 

1 / Q 



1 / 3 
2iVL V 2 1 



|a 3| 2 - 2a 30 ai2 + 2|a 12 | 2 ^ ^ |A| 4 A ds 
+ y* C( |A| 7 + (AiHEil + |Ex|) rfs. 

(5.29) 

Now define 

i = Ai - // 5 (Ai,0- (5.30) 
In terms of this new variable we have: 

A = A + / £*21 |i| 2 i + «32 |A| 4 i rfs 
•/ 

+ ^ ( ^|a 03 | 2 - 2a 30 a 12 + 2|a 12 | 2 ) |J |i| 4 i ds 

+ J* 0(\A\ 7 + |E| ) ds, (5.31) 
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where |E| is estimable in terms of |E| for \A\ < 1. This completes the proof. 

6. Large t behavior of solutions to the perturbed normal form equa- 
tions 

We now consider the large time behavior of solutions to the ordinary differential equations of the 
form (|5.2|) . In particular, we compare the behavior of solutions of fl5.2|) to those of the equation 
with E set equal to zero. 

Thus we consider the equations: 

P' = ic 2 i\P\ 2 P + (ic32- 7 ) \P\ 4 P + Q(*), (6.1) 
a' = ic2i\a\ 2 a + (ic 32 — 7) |a| 4 a, (6.2) 

where 7 > 0. We first consider the unperturbed equation, ( |6.2| ). Multiplication of ( p.2|) by a 
and taking the real part of the resulting equation yields the equation: 

r' = -2 7 r 3 , r = \a\ 2 . (6.3) 

Integration of (|6.3| ) yields: 

r 2 (t) = r 2 Q (l + 4 7 r 2 Q t)' 1 . (6.4) 

To prove the above lemma, we begin by multiplying ( |6.1|) by (3 and taking the real part of 
the resulting equation. This gives: 

r\t) = -2 7 r\t) + Q(t)£(t) + Q(t) (3(t), (6.5) 

which implies the differential inequality: 

r\t) < -2 7 r 3 (t) + 2 |Q(*)| r*(t). (6.6) 

We now prove 

Lemma 6.1. Suppose r(t) = \f3(t)\ 2 satisfies (\6.(\ ) with: 

|Q(t)| < Qo (t}- 5 ~ 5 , 5>0. (6.7) 

Then, 

< (l + 4|/3 |V) _1 (mi 4 + (6.8) 
where m* = max{l, 47|/? | 4 }- 
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proof: Use of ( |6.7| ) in flBTB] ) gives the inequality 

r\t) < -2 7 r 3 (t) + 2Q (t)~*~ S r*(t). (6.9) 
Multiplication by r(t) gives: 

< -47 z 2 (t) + AQ (t)-^ 8 z'i(t), where z(t) = r 2 (t). (6.10) 
Note that the equation 

C(t) = -4 7 C 2 W, C(0) = ^0 = |/5o| 4 (6.11) 

has solutions: 

C(t) = zo (1 + 4 W t)- 1 . (6.12) 
Anticipating this as the dominant behavior for large t, we define: 

z(t) = C(*) (6.13) 

Substitution into ( |6.10p and simplifying gives: 



5 r _3 



< : , ; U 7 ^(^-1) + C -^r (l+4zo7t)Ht)-~ d R*(t). (6.14) 
1 + 4 z >yt \zq\a 



We now consider the last term in ( |6.14[ ). Noting that 



(1 + t) 1 < m* (1 + 4-fz Q t) \ m* = max{ 1,472:0}, (6.15) 
we have for any e > 0: 



4 



, „ Q m$ emit) 

S t> i — r, r; x 



e |^ |i(t)f +5 (l + 4zo7*)* 

3 



gg ml £3 
S w — — 2 r - + ^2 



£ fz|(t} 1+ f 5 (l + 4z o7 t)- 

The last inequality follows from the inequality: ah < p~ 1 (ea) p + q~ l (b/e) q , p^ 1 + q^ 1 = 1, for 
I and 9=|- 



the choice p — | and g — 
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This last estimate can now be used in ( 6.14|) and implies: 



& < 



+ C 2 e* + 4 7 z 



< 



l + 4^o7t l + 4z ^t 

which when setting C%e* = 2720, is 

2720 R(R-2) + C Q ° m 



R(t) + Ci 



8 3 
, 5 



Qo ml 



£5 Z l ( t ) 



1+15 



i + 42 07 t ^ Zo (ty 

Now set i? = 2 + 5*. Since the term proportional to S 2 is negative, we obtain 



S' < 



1+4 z ^t 
Multiplication by 1 + 4 z^t yields: 



S + 



dt [ (1 + 4 z 7t) S(t) 



< C 



< c 



Qo 1 + 4 2 7t 1 



(1+*) (f)* 



7 5 ^o 



Integration of ( |6.18| ) from to t implies: 



S'(t) < C 



Qo ml 



7^0 

The Lemma now follows from (|6.19|) and the relation: 



=z{t) =c{t)R{t) = m (2 + s(t)). 



(6.16) 



(6.17) 



(6.18) 



(6.19) 



(6.20) 



7. Asymptotic behavior of solutions to the nonlinear Klein- Gordon 
equation 

In §2 we proved that local in time solutions, u(£), exist in the space C°(— T*,T*; X ), for some 
T*,T* > 0. We now study obtain the required 'a priori bounds to ensure (a) persistence of the 
solution, u(t), as a continuous function with values in Xo (T* = T* = 00) and (b) the decay 
of solutions as t — > ±00 in suitable norms. Due to the linear estimates of section 2, we require 



41 



more stringent hypotheses on uq and u\. These linear estimates require finiteness of W 1 ' 4 ^ 3 and 
W 1 ' 8 / 7 norms which, by interpolation, are controlled under the assumption Uo G X (see section 
2). Specifically, u G W 2 ' 2 n W 2 ' 1 and m G W 1 ' 2 n W 1 ' 1 . 

Using the results of the previous section, the original dynamical systems ( f4.1| ) can now be 
rewritten as: 

u(x,t) = a(t) tp(x) + r)(x,t), 
7)(x,t) = r)i(x,t) +T] 2 (x,t) +r) 3 (x, t), 
a{t) = A(t) e im + A{t) e~ im , 

A(t) = A(t)+h(A(t),t), (7.1) 

Here, h(A,t) is a smooth periodic function of t, cubic in A for A — > 0, rji and rj 2 are defined by 
( f4.12| )- ( |4.13| ), and A(t) satisfies the perturbed normal form equation: 

A' = i\c 21 \A\ 2 A - 7 |i| 4 i + i\ 2 c 32 \A\ 4 A + 0{\A\ 7 ) + E(t,a,rj). (7.2) 
Here E = E o (I + h), with E given by ( ^49|) . 

In Q, 

7 = 2 r = ^ (p c /,tf(B-3n)p c ^)>o. (7.3) 

The constants C21 and C32 are real numbers which are computable by the algorithm presented 
in section 5. 

The above definitions of A, A, and h, together with the estimates proved below, can be used 
to verify in a straightforward manner the assertions of Theorem 1.1 concerning R(t) = 
and 6(t) = arg A(t). 

To proceed with a study of the t — > 00 behavior, recall that: 

(d 2 + B 2 )r h = 0, m (x,0) = P c u , d tVl (x,0)= P cUl , (7.4) 

(d 2 + B 2 ) V2 = Aa 3 P c <^ 3 , V2 (x, 0) = 0, d tV2 (x, 0) = 0, (7.5) 
(d 2 + B 2 ) 773 = AP C (3aV 2 ^ + 3 mprj 2 + r? 3 ) , V3 (x, 0) = d m {x, 0) = 0. (7.6) 

To motivate the strategy, we first argue heuristically. Since 7 > 0, if E is negligible, then 
\A\ ~ (t) -1 / 4 as t -»• 00, and therefore by QOBD , a ~ (t) x / 4 , as t — > 00. It then follows from 
( |7.5| ) that, in appropriate norms, that rj 2 ~ (i)~ 3 / 4 and 773 ~ (t)~ 1+S for some 5 > 0. 
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To make all this precise requires d priori estimates on the system the above equations. We 



now recall Lemma 6.1 concerning equations of the form ( |7.2|) . This result gives conditions 
ensuring that A behaves like the solution of the equation obtained by setting E to zero. 

Our next task is to obtain an upper bound on E(t) in ( [4.49| ) of the form ( |6.7| ). Since the 
equation for A(t) is coupled to that of 77, the factor Q in ( |6.7D will depend on r\ and A. The 
next proposition, together with Proposition 4.7, provides estimates for the individual terms in 
E. 

It is convenient to introduce the notation: 

[A] a (T) = sup (f) a |il(f)| (7.7) 

0<t<T 

foU(T) = sup (t) a \\ V (t)\\v (7-8) 

0<t<T 

To avoid cumbersome notation, where it should cause no confusion, we shall abbreviate ex- 
pressions like [• • -](T) by [•••], until it is necessary to make the dependence on T explicit. In 
the estimates below, we shall often use the notation to denote a constant depending on 
some W k,p norm of the bound state, tp. Under our hypotheses, ip is a sufficiently smooth and 
exponentially decaying function for which these norms are finite; see [0. 

Because of Proposition 5.1, we need only estimate the terms of E, given in ( f4.49|) . The 
following proposition is the main step toward the estimate on E. 



Proposition 7.1. Estimates on the terms in E(£): 

For < s < t, the terms of E(s), as defined in ( |4.4fip , are estimated as follows: 



f m 2 


< 




f vrf 


< 





.,3/4 



] 8,3/4 
A "9/4 



-7/4 



(7.9) 
(7.10) 



(iii) \F 2 (a }Vl )\ < |A| [A\l /4 ||u ||x (t)- 13/8 , 

(iv) |F a (a, 773)1 < C V \X\ [A}j /4 C( [A] 1/4 , [rj] 8j3/4 , [B V3 U, 1/A+ao ,\\u \\x) (*)" 



(7.11) 



-5/4-CT0 



fvi) 



< 



C v |A| 2 \A \ 2 [A\l /4 (ty 



F 2 (a,C 2 ) < CJA| 2 [A}\ /4 [\A\M'h,A (t)~ 7/ \ 



vn 



rpnr 



< |A| C, {[W\A'\} 7/4 + [A]l /4 + [\A\\A 



|3\//| 



J 7/4 



-13/8 



(7.12) 
(7-13) 
(7.14) 

(7.15) 
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In (iv), C(ri, r 2 , r 3 , r 4 ) is bounded for Y^j \ r j\ bounded and tends to zero as J2j \ r j\ tends to zero; 
see Proposition 7.4. 

We now embark on the proof of this proposition. Parts (i) and (ii) follow by Holder's inequality. 
To prove part (iii), apply Holder's inequality and the linear propagator estimates of Theorem 
2.3. We now focus on (iv)-(vi). 

Estimation of F 2 (a,r) 3 ) 

Recall (see |4~l8p that 

F 2 (a,rj 3 ) = 3Xa 2 J <p 3 m . (7.16) 

We start by estimating H^Hs- We first express rj 3 , defined in (|7"Tj| ), as 

3 . 

r) 3 (t) = X / E^t - s)P c (3a 2 <f 2 n + 3a<^n 2 + r/ 3 ) ds, (7.17) 
3=1 Ij 

where 

/i = [0,t/2], I 2 =[t/2,t-l], and I 3 = [t-l,t]. (7.18) 
We estimate each integral separately. 

3 

\\V3(t)\\ 8 <\X\J2 / \\E l (t-s)P c (3a 2 ( p 2 r 1 + 3a^ 2 + r/ 3 ) || 8 ds. (7.19) 

7=1 Jl i 



1,8/7 ds, 

t-1 

l-fl N || , ,|l 

1.8/7 



The integrands are estimated as follows using the linear estimates of Corollary 2.1: 
/ \\E 1 (t-s)P c {---}\\ 8 ds < C f t/2 

I ii^(t- a )p c {...}ii 8 d8 < c r l \t- s \-^\\{---}\\ 1 ^ds, 

JI 2 Jt/2 

I \\E 1 (t-s)P c {-.-}\\ 8 ds < C f |t- s |- 3 / 8 ||{...}|| li8/7 ds, (7.20) 
Ji 3 Jt-i 

where 

ll{"-}||l,8/7< C(\A\ 2 \\<p 2 v\\l,8/7 + \A\ \\w 2 \\l,8/7 + \\V%,8/t) ■ (7.21) 



Proposition 7.2. 

(0 \A\ 2 \\<p 2 V \\ lj8/7 < C,\A\ 2 (\\n\ls + \\Brj\U), 

(ii) \A\ \\<prj%, 8 /7 < C V \A\ {\\n\\l + \\Brj\UM\,), 

(Hi) \\ V %, 8/7 < C(\\T,\hiM MIT (7-22) 
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proof: We prove part (i). The estimates (ii) and (iii) follows similarly. 

||<^||i,8/7 < C (\\(p 2 r]\\ 8/7 + \\(pd(prj\\ % / 7 + ||^ 2 ^|| 8/7 ) 

< C(l^ 2 ||4/3 IN Is + II^IU/3 IN Is + ll^l Is/5 H^IU), (7.23) 

by Holder's inequality. To express the right hand side of (|7.23 ) in terms of ||5r/||4, and thereby 
completing the proof of part (i), it suffices to show that: 

Ml* < C\\B V \\ A . (7.24) 

This follows if we show that the operator d^B" 1 is bounded on LP (with p — 4). 

To prove the LP boundedness of d^B^ 1 , we can apply the results on the wave operator, W* 
in §2.2. Indeed, for any g G L p , we have using the boundedness of the wave operators on W 1,p , 
for p > 1, that 

WdiB-^Wp = WdiW+B^WZgWp 

< C\\W + B, l Wlg\\ w ^ 

< C\\B^Wlg\\ w ^ 

< C\\Wlg\\ p < C\\g\\ p . 

Remark: We offer here an alternative proof which does not make use of the wave operators, 
and therefore applies under weaker hypotheses on V. We begin with the square root formula 
H and (p|): 



sec 



For any i/j G V(B 2 ): 

B- 1 = 7T- 1 / w- 1/2 {B 2 + W)- 1 dw. (7.25) 

J 

By (|7.25 ) and the second resolvent formula we have: 

poo 

diB- 1 = d^ 1 + it- 1 / w- 1/2 dABl+w)- 1 V (B 2 + w)- 1 dw. (7.26) 

Jo 

The boundedness in LP (p > 1) of the operator diB^ 1 holds because £i(|£| 2 + m 2 ) -1 / 2 is a 
multiplier on L p ; see HS4| . Similarly, di(B 2 + w)' 1 is bounded on L p for any w > and 
estimation of the second term in (|7. 26|) is reduced to estimation of the norm of the operator 
V (B 2 + w)- 1 . Note that: 

roo „ 

(B 2 + w)- 1 = / e' t{B +w) dt (7.27) 
Jo 



Since inf a(B 2 ) = fi 2 > 

n „-t(B 2 -n 2 /2)\, 

\B{LP) 



e -* (B 2-n72)|| < c (7 _ 28 ) 
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and therefore 

poo 

\\{B 2 + w)- l \\ B[LP) < C J e ~ t{n2/2+w) < C\l+w)-\ 
Boundedness in L p of diB~ x now follows. Namely, 



(7.29) 



l B- 1 \\ Bm < + C'WdiiBt + wy^BwWVWn / w- l '\l + w)- 1 dw < 00. 

** 

(7.30) 



From ( |7.22| ) we see that H-B77H4 must be estimated. Recall that 

Br] = B% + Bt] 2 + Br] 3 . 

By@ 



and so 



Brn = E {t)BP c u + E 1 (t)BP c u 1 , 

\\B Vl (t)\U < C(t}-^ 2 Huollx; 

see ( p.7| ) for definition of ||u || X - 

By© 

B V2 = A / EAt - s) a 3 (s) BP c p 3 ds 
Jo 

which can by Theorem 2.3 be estimated as 

\\B V2 (t)\U < C\X\ f\t-s\-^ \A{s)\* ||P c ^ 3 || 4/3 ^ < C v \\\ [Af 1/4 (t)- 1 2 

J 

The estimation of [ | ^773^) | (4 is more involved. By (|7 



Br] 3 = A / sin B(t - s) P c 3<r ip 2 77 + 3a iprf + rf 
Jo L ' 



ds. 



By Theorem 2.3, 



(7.31) 

(7.32) 
(7.33) 

(7.34) 
(7.35) 

(7.36) 



\\Bvs(t)\U < C\X\ / \t - s^ 2 \ \A\ 2 ||AIIm/3 + \A\ iWWws + Wv'Wia/s] ds. (7.37) 



The following proposition provides estimates for the integrand. It is proved in the same manner 
as Proposition 7.2. 

Proposition 7.3. 

\A\ 2 \\<p 2 v\\i A /3 < C v \A\ 2 (\\ V \\s + \\B Vl \U + \\B V2 \U + \\B V3 \U), 



(7.38) 



\A\ ||w 2 ||i,4/3 < C v \A\ (Ml + Mis \\Br,\\< 

\\v%A/s < 3 iMil ||^|| 2 . 
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Anticipating the behavior 

\A(t)\~r l /\ |to(t)|| 8 ~ r 3/4 , 

and using ( |7.33| ), we find that Bf] 3 is driven by terms which are formally of order (t) -1 . Esti- 
mation in L 4 will lead to convolution of (t)^ 1 with (t)^ 1 ^ 2 giving a rate of (t)~ 1//2+<5 °, for any 
5 >0. 

For < t < T, with T fixed and arbitrary, we have: 



l^l 2 11^114/3,1 < C, [A) 2 1/4 (W 8i3 /4(t)- 5/4 + [B^y.it)- 1 

+ [BmkiMt}- 1 + [B V3 \± jl/2 „ So (t)- 1+ ^ 



where 5 is positive and arbitrary Also, 



\A\ \\vv 2 \\ia/3 < C v [A] 1/A ( [rf\l m {t)~ 7 " + feN^i)" 1 

11^111,4/3 < lhlil.2 W| 3 /4 (*)~ 3/2 - 



It follows that for < t < T 



HS773II4 < C„|A| f \t - si' 1 ' 2 ds(t)- 1+So G(A, V ,T), (7.39) 
Jo 

where 

G{A,r],T) = [A] 2 l/4 (Ms.3/4 + [^1)4,1/2 + [5r/ 2 ] 4 ,i /2 + [fltfeki/a-*, 

+ [4 1/4 ( Ws.3/4 + M 8,3/4 [BVlUfi ) + II^Hl,2 Ws,3/4- 

Therefore, for < t < T: 

\\Bm{t)\U < ^\X\ G(A, V} T) (ty 1 ^ 50 , (7.40) 

for any 5 > 0. We now use the estimate ( |7.4(J| ) in ( |7.22| ) in order to obtain a bound on | \r] 3 (t) | | 8 . 
First, a consequence of ( |7722| ), QT53D , ( |7^5Q and (|7^D|) is that for < t < T: 



|^| 2 ||^Hl,8/7 < ^ [A] 2 /4 ( M 8i3/4 + Huollx + [AW,, + [B V3 ] 4>1/2 ^ ) (t) 

\A\ ||^ 2 ||l,8/7 < C v [A] 1/4 ( [< 3 /4 + {||u ||x + [S77 3 ]4,l/2-«5o + [4l/J M 8,3/4 ) (*> 

ll^lks/r < CV(|M|i j2 ) MjJ /4 (*)- 5/4 . (7.41) 



l+5 
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Substitution of (fT41|) into ( fTT7| - [r20| ) leads to an estimate for H^^Hs- For <t <T: 



\\m(t)\\s < \\\ (t}- 1+s » ( [A]? /4 {[77] 8j3 /4 + ||u ||x+ [A]l /A + [BmU >1/2 „ So } 

+ Hl/4 { Ms,3/4 + [ Kl|x+ [5t7 3 ] 4,l/2-<5 

+ G(iNii,2^) ), (7-42) 

where #o is arbitary. 

Finally we can now estimate ^2(0,773) using ( [4.18| ) and ( |7.42| ). Choosing 5q so that 



we have 



— 1 + 5 = —3/4 — cr , with cr > 0, 



Proposition 7.4. 

1^2(0, 773)1 < C V \X\ [A] 2 1/4 ( [A]~y 4 { [77] B ,3/4 + ||U |X + [A]? /4 + [5^4,1/4+ao } 

+ [^]l/4 { Ws,3/4 + [ll u o||x + [A]l /4 + [ J Br7 3 ]4,l/4+a ] foks/J 

+ C{\\r,\\w,<p) [ V }t% 4 ) (t)" 5 / 4 --. (7.43) 



Estimation of F 2 (a,r]™ rl ) and F 2 (a, 77" 



Jir2\ 



By (pi 



A 
2 

3A 



- 3fi + iO)" 1 e*"* P c ^, and 



q \ ft 

: _^ 5-1(5 - 3ft + 70)- 1 / e ii? ('- s ) e i3Cs A 2 A' dsP, 

2 jo 

Estimation using Proposition 2.2 gives: 



3 |A| \a\ 2 
C\X\ 2 \A \*\A\ 2 



- 1 e wt P c ^ 



< C\\\ 2 \A \ Z \A\ 2 || (x)V || 2 H^-^S-Sfi + ^O)- 1 e wt P c fl-y|| 2 

< CJA| 2 |A | 2 [A]\ /A (t)-l+ & . (7.44) 

Here, we have used that 1/2 + 6/5 = 5/4 + 5, for some 5 > 0. The constant depends 
on || (x) cr B~ 1 (x)~ cr \\B(L 2 )- It is easy to check that this norm is bounded if we replace B^ 1 by 
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B 2 . The estimate of interest is reduced to this case using the Kato square root formula (|7.25|) . 
Similarly, we have 



F 2 (a,T].. 



= 3|A| \a\ 2 
< C|A| 2 U| 2 



^T 2 



J if 3 [B(B - 3ft + is)]- 1 J e i(t-»)(B-xi) A 2 A , ds p 

< C^\\\ 2 \A\ 2 (\t-s)-$+» \A\ 2 \A'\ ds 
Jo 



< C v \X\ 2 [A]U [\A\ 2 \A'\] 5/4 



(7.45) 



where we have taken // sufficiently small so that the last integral t— integral is convolution with 
an L 1 function, which then preserves the decay rate , (t)~i, which exceeds (t)~^ 

Estimation of E%J : 



To estimate (|Q7D we use Propositions 2.1 and 2.2. Due to the singularity in the resolvent 
at frequency 3fl, the case j = 7, is most difficult and we focus on it. To treat this singularity, 
we use Proposition 2.2, for n = 3: 



\\{ x )-°(B-C 7 )- 2 e m - s) {x)-°n2 < C(t-sY 
Use of this estimate in ( |4.47| ) yields 



1,2- 



(7.46) 



\E%\ < 

C V \X\ 2 ( [1^1 V'l]r/4(r 7/4 + [A] 2 1/4 (tr 13/8 + U\ 4 \A"\ + \A\ 3 \A>\ 2 } 9/ ,(t)- 13 / 8 ) , 

from which estimate (vii) of Proposition 7.1 follows. This finally completes the proof of Propo- 
sition 7.1. 

Proposition 7.1 can now be used to obtain the desired estimate for E(t) and therefore A(t), 
for \A\ sufficiently small. Note that the right hand side of the estimates depend on \A'\ and \ A"\. 
These can each be estimated directly from the equation for A, ( 4.48|) , and its derivative. The 
result is that \A'\ and \A"\ may be estimated by [A]w 4 (i)~f , the only effect being a change in 
the multiplicative constants which is independent of A and rj. This observation together with 
Proposition 7.1 yields: 

Proposition 7.5. There is positive number 8 such that 



\E(t)\ < Q (A,rj)(t)-r 



(7.47) 
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where 

Qo{A,r]) = [77] 1,3/4 + [4 1/4 (l + C ([A]i /4 , M8.3/4, [Br] 3 U ;l/4+ao , ||u ||x ) ) (7.48) 
It now remains to estimate |H|s- 



IMIs < IMIs + \\v2Ws + WmWs 

< \\E (t)P c u \\ 8 + ||Si(t)PcMi|| 8 + |A| f t \\E 1 (t-s)a 3 (s)P c <p% ds + \\ V3 \\ 8 . 

Jo 

Using Theorem 2.3 and the bound (|7.42|) we get 
Proposition 7.6. 



< C v ( ||u ||x + [A] 3 1/A 

+ [Al/d M 8,3/4 + Kllx + [A]l /4 + [£r/3k 1/4+ff0 } 
+ [A}i/4 { [v]l,3/4 + lluollx + [A]? /4 + [B V3 ] 

4,1/4+c 

+ c-(|M| lj2 ^) Mgf /4 ) <r 3/4 - 



W 8,3/4} 



(7.49) 



The following proposition summarizes our labors. 
Proposition 7.7. For any T > 0: 



Ws,3/4(T) 


< 


^(||u ||x + [A] 3 1/4 (T) + 


[^12,1/4^ (r) 




+ 


cA\\v\\i,2) [v]%Ut)) 






< 


(\A \* + Qo(A, V )i ) 




[^ 3 ]4,l/4+ CT0 ( T ) 


< 


(M1/4CO + + 


P^lli/aCO + [S7 73 ]l, 1 /4 +CT0 (T) ) 


SUP \\v(t)\\l,2 


< 


C £(w ,«i) < C||m ,mi||i, 2 . 





0<t<T 

The first three estimates are proved above while the last follows from conservation of energy 
(see section 1) and the decomposition of the solution. 
Now define 

M(T) = W 8 , 3/4 (T) + [A] 1/4 (T) + [B V3 U il/4+l70 (T). (7.50) 
Then, combining the estimates of the previous proposition we have, for some a > 0: 



M(T)( 1-M(T) Q ) < CL ||u Q | 



(7.51) 
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If M(0) and ||uo||x are sufficiently small, we have by the continuity of M(T), that there is 
a constant M*, which is independent of T, such that for all T 

M(T) < M*. (7.52) 

This completes the proof of Theorem 1.1. 



8. Summary and discussion 



We have considered a class of nonlinear Klein-Gordon equations, (1.1), which are perturbations 



of a linear dispersive equation which has a time periodic and spatially localized (bound state) 
solution. The unperturbed and perturbed dynamical systems are Hamiltonian. We have shown 
that if a nonlinear resonance condition ( |1.8|) holds, then solutions with sufficiently small initial 
data tend to zero as t — > ±oo. This resonance condition is a nonlinear variant of the Fermi golden 
rule (|1.22j ). A consequence of our result is that time-periodic and spatially localized solutions do 



not persist under small Hamiltonian perturbations. Time-decay of small amplitude solutions is 
also a property of the translation invariant (V = 0) nonlinear Klein-Gordon equation. However, 
the presence of a bound state of the unperturbed problem, causes a nonlinear resonance leading 
to the anomalously slow radiative decay of solutions. 

We now conclude with some further remarks on the results of this paper, mention directions 
currently under investigation and some open problems. 

1. Anomalously slow time-decay rates: It is natural to compare the time-decay rate of 
solutions described by Theorem 1.1 with those of related problems. 

(la) Translation invariant linear Klein-Gordon equation (V = and A = 0): 
We shall refer to free dispersive rates of decay as those associated with the constant coefficient 
equation: 

d 2 t u - Au + m 2 u = 0. (8.1) 

Results on this are presented in §2. Roughly speaking, if the initial data has a sufficient number 
of derivatives in L p , 1 < p < 2, then the solution decays at a rate 0{t~ n ^~~i> r ' > ) in LP' . Here, 
p- 1 + (p'Y 1 = 1. 

(lb) Translation invariant nonlinear Klein Gordon equation (V = and A ^ 0): 
For small initial conditions it has been shown that solutions decay at free dispersive rates; see 



65 and references cited therein. 
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(lc) Linear Klein-Gordon equation with a potential having a bound state, as hypothesized 
(1^/0 and A = 0): By the spectral theorem, a typical solution will decompose into a linear 

superposition of (i) a bound state part, of the form: R cos(Qt + p )ip(x), and (ii) a part which 
disperses to zero at free dispersive rates. 

(Id) Nonlinear Klein-Gordon equation with a potential having a bound state, as hypothesized 
(V ^ and A ^ 0): Whereas the decaying part of the solution in the above cases tends to 
zero at a free dispersive rate, Theorem 1.1 implies that the decay rate is anomalously slow. In 
particular, the decay rate obtained in L 8 is 0(t~^), while the free dispersive rate in L s is 0(t~s). 
This slow rate of decay, due to the nonlinear resonant interactions give rise to a long-lived or 
metastable states. 

2. Dispersive Hamiltonian normal form 

In a finite dimensional Hamiltonian system of one degree of freedom, the general normal 
form is: 

A' = i (do + c 21 \A? + cs2\A\ 4 + ... + c n+1 , n \A\ 2n + ...) A, (8.2) 

where c nj „ + i are real numbers. In the present context we have the dispersive Hamiltonian normal 
form 

A' = (k 10 + k 21 \A\ 2 + hMf + ... + k n+1<n \A\ 2n + ...) A, (8.3) 

where 

kn+l,n ^n+l,n ^Cn+l,ra 

are, in general, numbers with real and imaginary part. Contributions to the real parts of 
coefficients come from resonances with the continuous spectrum. If ( p..8|) fails there are two 
possibilities; either (a) 3fl does not lie in the continuous spectrum of B (it lies in the gap 
(0, m)), or (b) 3f2 lies in the continuous spectrum of B but we are in the non-generic situation 
where the projection of ip 3 onto the generalized (continuum) eigenmode of frequency 3f2 is zero. 
In either case, we expect that, typically, internal dissipation would arise in the normal form at 
higher order. More precisely, we conjecture that the leading order nonzero d nt+ i t n tl which would 
correspond to a resonance with the continuum of a higher harmonic: G a con t(B), > 3 
is always negative. (If the sign of d nt+ i trit were positive, this would seem to be in violation 
of the conservation of energy and the implied Lyapunov stability of the zero solution.) The 
corresponding decay rate would then be slower, specifically 0{\t\~^). 

From this perspective one may expect the existence of breather solutions of integrable non- 
linear flows like the sine-Gordon equation or the extremely long-lived breather -like states of 
the 4 - model as corresponding to the case where to all orders the coefficients d n+ i >n are zero. 
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Does the vanishing of all such coefficients have an interpretation in terms of the infinitely many 
time-invariants for the integrable flow? 

3. Multiple bound state problems: Of interest are problems where the underlying potential, 
V(x), supports more than one bound state. What is the large time behavior of such systems? 
Our analysis and the above remarks suggest that a more general normal form could be developed, 
and the corresponding ( in general slower) decay rates anticipated. Is it possible that a resonance 
among the multiple "discrete oscillators" can be arranged so that one gets a persistence of 
nondecaying solutions, or does radiation to the continuum always win out? 

4. Relation to center manifold theory: Our program of decomposing the original conserva- 
tive dynamical system into a finite dimensional dynamical system ( 1.20Q which is weakly coupled 
to an infinite dimensional dynamical system is one commonly used in dissipative dynamical sys- 



tems. There, it is often possible using the center manifold approach |14 |, [ pTjl , to construct an 



invariant center-stable manifold. The dynamics in a neighborhood of an equilibrium point are 
characterized by an exponentially fast contraction on to the stable-center manifold. The con- 
traction is exponential because the part of the linearized spectrum associated with the infinite 
dimensional part of the dynamics is contained in the left half plane. In the current context of 
conservative dynamical systems, the linearization about the equilibrium point (here u = 0) lies 
on the imaginary axis; in particular, two complex conjugate eigenvalues and continous spec- 
trum from ±im to ±ioo. The analogue of dissipation is the mechanism of dispersive radiation 
of energy, related to the continuous spectrum, and the associated algebraic time-decay. 

There is recent work on the application of center manifold methods to certain special conser- 
vative dynamical systems of nonlinear Schrodinger type; see the center manifold analysis of 



applied to problem studied by the authors in |JIJ. It would be of interest to understand whether 
geometric insight on the structure of the phase space for problems of the type considered in this 
paper can be obtained using the ideas of invariant manifold theory. 

5. Systems with disorder: In this paper, we have seen the effect of a single localized 
defect on the wave propagation dynamics in a nonlinear system. Of great interest would be an 
understanding of the effects of a spatially random distribution of defects modeled, for example, 
by random potential V(x) on the localization of energy in nonlinear systems such as (|1 . 1|) . 
Related questions are studied in |p2|, f25fl, [TO . 
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